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I.  INTRODUCTION 


This  final  report  covers  the  three  year  contract  period  1  January 
1977  through  December  31,  1979.  During  this  time  research  work  has  been 
carried  out  by  Professor  A.E.  Siegman  in  the  areas  of  mercury  lasers, 
picosecond  spectroscopy,  and  optical  beam  and  resonator  calculations; 
and  by  Professor  R.L.  Byer  in  the  areas  of  optically  pumped  lasers; 
two-photon  laser  systems,  optical  parametric  amplifier  oscillator  studies, 
and  stimulated  rotational  Raman  scattering.  During  the  past  three  year 
contract  period  important  results  have  been  obtained,  including: 

A.  Optically  Pumped  Atomic  Mercury  Laser 

B.  Picosecond  Spectroscopy  Experiments 

C.  Fast  Transform  Calculation  Methods  for  Optical  Beam 
and  Resonator  Calculations 

D.  Calculation  Methods  for  Nonlinear  Wave  Interactions 

E.  Basic  Laser  Resonator  Theory 

F.  Diffraction-Biased  Ring  Unstable  Resonators 
Phase  Conjugate  Mirror  Resonators 
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J.  Frequency  Conversion  via  Four  Wave  Mixing 

K.  Optical  Parametric  Oscillator  and  Optical  Parametric 
Amplifier  Studies 

L.  Efficient  Frequency  Conversion  via  Stimulated  Raman  Scattering. 

M.  Current  Research. 
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Details  of  these  results  are  given  below  together  with  the 
resulting  publications  and  the  five  Ph.D.  theses  which  have  been  com¬ 
pleted  during  this  period. 

II.  REVIEW  OF  ACCOMPLISHMENTS 

A.  Optically  Pumped  Atomic  Mercury  Laser 

During  the  early  part  of  this  program  we  were  engaged  in  developing 
under  AFOSR  support  a  new  optically  pumped  mercury  laser  oscillating  at 

O 

5461  A  in  the  green.  This  laser  oscillates  in  a  laser  tube  containing 
a  mixture  of  nitrogen  and  mercury,  when  pumped  by  a  special  high-bright¬ 
ness  mercury  discharge  lamp.  Because  this  laser  involves  no  noisy  plasma 
discharge  in  the  laser  tube  itself,  and  because  it  oscillates  at  a  con¬ 
venient  green  wavelength,  it  is  of  particular  interest  as  a  possible 
laser  gyro  gain  medium;  as  a  visible  wavelength  standard;  and  possibly  as 
a  general  purpose  low-power  green  laser. 

During  the  first  year  of  the  program  reported  here,  the  work  on 
this  project  was  brought  to  a  successful  conclusion  and  then  terminated. 

A  special  diffusion-bypass  low-pressure  mercury  arc  lamp  to  meet  the 
pumping  needs  of  this  laser  with  minimum  mercury  fill  requirements,  suit¬ 
able  for  single-isotope  operation,  was  invented,  constructed,  and  operated 
very  successfully  for  many  hundreds  of  hours.  Parametric  studies  of  the 
laser  itself,  including  measurements  of  gain  and  power  output  as  functions 
of  pressure,  bore  diameter  and  pumping  intensity,  were  carried  out. 

Because  this  laser  is  of  particular  interest  as  a  visible  wavelength 
standard,  measurements  of  saturated  absorption  in  iodine  using  this  laser 
formed  one  of  the  major  goals  of  this  program.  These  measurements  were 
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carried  out,  and  a  strong.  Isolated  saturated-absorption  feature  In 

127 

I2  was  observed.  The  pressure  shift  of  the  laser  transition  against 
this  iodine  resonance  was  measured,  and  it  was  found  possible  to  center 
the  Hg  laser  transition  exactly  on  the  iodine  resonance  by  varying  the 
nitrogen  pressure  in  the  laser  tube.  These  results  laid  a  firm  founda¬ 
tion  for  possible  later  absolute  wavelength  stabilization  of  this  laser 
against  a  low-pressure  iodine  cell. 

Work  on  this  topic  under  AFOSR  support  was  reported  in  detail  in 
the  following  publications: 

-  Max  Artusy  and  A.  E.  Siegman,  "Saturable  Absorber  Overlap  of 
Iodine  127  with  the  Optically  Pumped  546  nm  Mercury  Laser," 
Appl.  Phys.  Letters  31,  333-334  (1  September  1977). 

-  Max  Artusy,  "Optically  Pumped  Mercury  Laser  for  Absolute 
Frequency  Stabilization,"  Ph.D.  Dissertation,  Department  of 
Electrical  Engineering,  Stanford  University,  Stanford, 
California,  1977. 

B.  Picosecond  Spectroscopy  Experiments 

Some  very  interesting  measurements  of  ultrafast  relaxation  times 
in  solids  were  completed  during  this  program,  using  the  novel  transient 
laser- induced  grating  technique  for  picosecond  spectroscopy  developed 
in  our  group.  The  measurements  demonstrated  in  particular  a  novel  photo 
acoustic  pressure- induced  amplitude  grating  effect  in  molecular  crystals 
of  a  type  which  does  not  appear  to  have  ever  been  observed  previously. 
This  work  was  reported  at  the  1978  OSA  Meeting  and  in  a  journal  article, 
of  which  a  reprint  is  attached  as  Appendix  I  of  this  report,  as  well  as 
in  a  Ph.D.  dissertation: 
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-  J.  R.  Salcedo  and  A.  E.  S legman,  "Laser  Induced  Photoacoustic 
Grating  Effects  In  Molecular  Crystals,"  IEEE  J.  Quant.  Electr. 
QE-15,  250-256  (April  1979). 

-  J.  R.  Salcedo,  "The  Picosecond  Transient  Grating  Technique: 
Studies  on  Energy  Transport  and  Photo-acoustic  Effects," 

Ph.D.  Dissertation,  Stanford  University,  Stanford,  California 
(February  1979) . 

-  J.  R.  Salcedo  and  A.  E.  Siegman,  "Amplitude  Grating  Effects 
Produced  by  Picosecond  Laser  Produced  Acoustic  Waves,"  (Paper 
THC10,  abstract  only),  J.  Opt.  Soc.  Am.  68,  1412  (October  1978). 
Paper  presented  at  the  1978  Annual  Meeting  of  OSA. 

Work  on  this  type  of  picosecond  measurement  using  the  transient  grating 
approach  is  continuing  in  Professor  Siegman* s  laboratory,  and  has  been 
adopted  by  other  groups  elsewhere. 

C.  Fast  Transform  Calculation  Methods  for  Optical  Beam  and 
Resonator  Calculations 

One  of  our  primary  research  areas  has  been  the  development  and 
verification  of  efficient  and  accurate  numerical  methods  for  handling 
optical  beam  and  resonator  calculations  that  are  amendable  to  purely 
analytic  treatment.  This  includes  many  situations  of  importance  in 
high  power  and  high  energy  lasers.  One  of  our  most  important  contri¬ 
butions  has  been  a  quasi  Fast  Hankel  Transform  technique  for  handling 
beam  calculations  in  cylindrical  coordinates  (in  contrast  to  Fast  Fourier 
Transform  techniques,  which  apply  in  rectangular  coordinates). 

In  addition  to  the  first  publication  of  this  transform  algorithm, 
during  the  contract  period  substantial  steps  were  accomplished  in 
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completing  and  testing  the  novel  Fast  Hankel  transform  algorithm  that 
we  first  invented,  and  in  applying  it  to  optical  calculations,  particu¬ 
larly  to  beam  and  resonator  calculations.  This  new  FHT  algorithm, 
first  published  under  AFOSR  support  in  mid-1977,  has  since  aroused 
wide  interest  and  been  used  in  several  other  AFOSR-supported  industrial 
high-power  laser  programs. 

Improvements  to  the  original  FHT  algorithm  under  this  program 
during  the  past  year  have  included  a  simple  "lower  end  correction" 
technique  to  compensate  for  the  finite  truncation  of  the  transform  integral 
at  the  lower  end  of  the  integration  range;  more  accurate  ways  of  handling 
the  discrete  Fast  Fourier  transform  steps  involved  in  the  calculation; 
and  better  understanding  of  the  choice  of  sampling  Intervals  in  the 
calculation.  These  improvements  brought  the  accuracy  of  the  transform 
close  to  the  limits  of  single  precision  arithmetic  on  our  machine  (frac¬ 
tional  error  ~  10  b)  .  In  addition,  a  two-stage  FHT  algorithm  requiring 
roughly  twice  the  storage  and  computation  time  was  developed  which  reduces 

the  error  of  the  algorithm  itself,  within  its  range  of  validity,  to  below 
-12 

~  10  which  is  probably  better  than  will  ever  be  needed. 

Publications  that  have  appeared,  or  are  in  the  final  stages  of 
preparation,  covering  the  FHT  algorithm  include: 

-  A.  E.  Siegman,  "Quasi  Fast  Hankel  Transform,"  Optics  Letters  JL, 
13-15  (July  1977). 

-  Shinan-Chur  Sheng,  "Studies  of  Laser  Resonators  and  Beam  Propa¬ 
gation  Using  Fast  Transform  Methods,"  Ph.D.  Dissertation,  Depart¬ 
ment  of  Applied  Physics,  Stanford  University,  Stanford,  California 
(1980). 
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A.  E.  Siegman  and  Shinan-Chur  Sheng,  "The  Fast  Hankel 
Transform  Algorithm,"  manuscript  in  preparation  0-980). 


D.  Calculation  Methods  for  Nonlinear  Wave  Interactions 

Another  result  of  our  program  has  been  the  development  of  an 
efficient  procedure  for  analyzing  important  nonlinear  interactions,  such 
as  harmonic  generation  and  parametric  or  Raman  amplification,  including 
realistic  considerations  such  as  arbitrary  beam  profiles,  diffraction 
and  depletion  effects,  phase  mismatch  and  walk-off,  though  not  (for  the 
minute)  time-varying  effects.  An  efficient  procedure  for  analyzing 
such  phenomena  is  important  because  these  interactions  provide  an 
important  way  in  which  the  energy  of  an  efficient  laser  beam  at  one  wave¬ 
length  can  be  converted  into  useful  and  powerful  radiation  at  other  wave¬ 
lengths,  either  longer  or  shorter,  where  an  efficient  laser  system  may 
not  be  available.  One  immediately  important  example  is  harmonic  genera¬ 
tion  and  parametric  amplification  with  strong  depletion,  using  an  unstable 
reasonator  beam  as  the  primary  source,  since  this  has  already  become  an 
important  practical  technique  for  generating  new  coherent  wavelengths. 

We  have  therefore  developed  and  made  preliminary  use  of  an  efficient 
numerical  method  that  accomplishes  these  objectives.  To  carry  out  the 
nonlinear  interaction  calculations  we  Fourier  or  Hankel  transform  the 
relevant  nonlinear  wave  equation  in  its  transverse  coordinates,  and  obtain 
coupled  differential  equations  of  the  form  dE^/dz  =  F^[E^]  ,  where 

=>  E^s.z)  is  a  spatially  transformed  wave  at  the  i-th  frequency  (i>^ 
and  s  is  the  transverse  transform  corrdinate(s) .  Optical  diffraction 
effects  can  be  entirely  transformed  away  in  this  transformation,  so  that 


the  function  F^[E^]  contains  only  the  nonlinear  interaction  of  the  wave 
with  all  other  components.  Fast  Fourier  or  Hankel  transform 

methods  are  used  to  evaluate  the  derivative  functions  F^[E^]  »  an<* 
standard  predictor-corrector  ordinary  differential  equation  (ODE)  methods 
are  used  to  step  foreward  along  z  . 

A  preliminary  account  of  this  work  was  presented  in  1977,  and  a  full 
published  report  is  now  in  press: 

-  A.  E.  Siegman  and  Shinan-Chur  Sheng,  "Optical  Beam  Calculations 
Using  the  Fast  Hankel  Transform,"  (Paper  FD10,  Abstract  only), 

J.  Opt.  Soc.  Am.  j>8,  1436  (October  1978).  Paper  presented  at 
the  1978  Annual  Meeting  of  the  OSA. 

-  Shinan-Chur  Sheng  and  A.  E.  Siegman,  "Nonlinear-Optical  Calcu¬ 
lations  Using  Fast-Transform  Methods:  Second-Harmonic  Generation 
with  Depletion  and  Diffraction,"  to  be  published  in  Phys.  Rev.  A 
(1980) . 

E.  Basic  Laser  Resonator  Theory 

Despite  the  now  rather  advanced  state  of  laser  theory,  several  very 
fundamental  aspects  of  laser  resonators  have  until  recently  remained 
unclear  or  uncertain.  For  example,  because  the  fundamental  propagation 
kernel  in  optical  resonator  theory  is  not  a  hermitian  operator,  the 
eigenmodes  of  optical  resonators  are  not  normal  (i.e.,  orthogonal)  modes, 
and  the  orthogonality  properties  of  optical  resonator  modes  have  remained 
uncertain.  All  laser  analyses  to  date  have  nonetheless  obtained  laser 
cavity  equations  of  motion  by  assuming  a  hypothetical  set  of  lossless 
and  orthogonal  cavity  modes,  despite  the  fact  that  real  optical  resonator 
modes  have  significant  diffraction  losses  and  are  not  orthogonal. 
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All  of  these  fundamental  questions  have  been  essentially  resolved 
under  this  program  during  the  past  year.  The  orthogonality  properties 
of  optical  resonators  have  been  found  to  be,  in  essence,  that  each  eigen¬ 
mode  going  in  one  direction  in  an  optical  cavity  is  orthogonal  to  all  the 
other  eigenmodes  of  an  adjoint  (biorthogonal)  set  of  modes  going  in  the 
opposite  direction  in  the  same  cavity.  In  addition,  a  rigorous  forip  of 
laser  analysis  has  been  developed  which  gives  laser  cavity  equations  of 
motion  in  terms  of  the  real,  lossy,  nonorthogonal  modes  of  the  real  laser 
cavity.  This  new  analysis  reduces  to  previous  results  in  the  appropriate 
limits,  but  is  much  better  adapted  to  handle  such  topics  as  high-gain 
lasers  with  rapidly  changing  mode  properties;  lasers  with  external  signal 
injection;  and  lasers  with  large  spontaneous  emission  noise. 

All  of  these  results  are  reported  in  detail  in  the  publications: 

-  A.  E.  Siegman,  "Orthogonality  Properties  of  Optical  Resonator 
Eigenmodes,"  Optics  Commun.  31,  369  (1979). 

-  A.  E.  Siegman,  "Exact  Cavity  Equations  for  Lasers  with  Large 
Output  Coupling,"  to  be  published  in  Appl.  Phys.  Letters  ,36 
(1980) . 

-  A.  E.  Siegman,  "Cavity  Equations  for  Laser  Cavities  with  Large 
Output  Coupling,  Injected  Signals,  and  Spontaneous  Emission 
Noise,"  manuscript  in  preparation  (1980). 

F.  Diffraction-Biased  Ring  Unstable  Resonators 

A  novel  type  of  diffraction-biased  ring  unstable  resonator,  with 
potential  application  to  ring  laser  gyros,  was  invented  and  analyzed 
under  this  program.  This  result  was  published  in: 


-  Shinan-Chur  Sheng,  "Diffraction-Biased  Unstable  Ring  Resonators 
with  Possible  Applications  in  Laser  Gyroscopes,"  IEEE  J.  Quant. 
Electr.  QE-15.  922-926  (September  1979). 

G.  Phase  Conjugate  Mirror  Resonators 

One  of  the  most  interesting  developments  in  optics  in  recent  years 
has  been  the  emergence  of  the  phase  conjugate  mirror,  using  four-wave 
mixing,  as  a  new  fundamental  optical  component.  During  the  past  year, 
we  have  carried  out  a  basic  analysis  of  the  resonator  mode  properties 
of  optical  cavities  having  a  phase  conjugate  mirror  at  one  end.  These 
results  are  reported  in  detail  in: 

-  Pierre  A.  Belanger,  Amos  Hardy,  and  A.  E.  Siegman,  "Resonant  Modes 
of  Optical  Cavities  with  Phase-Conjugate  Mirrors,"  to  be  published 
in  Appl.  Optics  19^  (1980). 

-  Pierre  A.  Belanger,  Amos  Hardy,  and  A.  E.  Siegman,  "Resonant 
Modes  of  Optical  Cavities  with  Phase  Conjugate  Mirrors:  Higher- 
Order  Modes,"  to  be  published  in  Appl.  Optics  19  (1980). 

H.  Hg^  Optical  Pumping  Studies 

Research  in  optically  pumped  lasers  occupied  most  of  our  effort  in 

the  first  year  of  this  program.  The  work  built  on  our  earlier  success 

with  the  demonstration  of  the  first  electronic  transition  optically 

pumped  laser  in  molecular  ^  We  followed  the  I.  work  by  demonstrat- 

2 

ing  an  optically  pumped  Na2  molecular  laser.  Those  early  experiments 
did  much  to  open  the  field  of  optically  pumped  lasers  such  that  operation 
of  a  wide  variety  of  lasers  of  this  type  has  now  been  accomplished. 

Under  this  program  we  extended  our  studies  to  Hg£  because  of 
its  potential  as  an  energy  storage  laser  with  tunable  output.  Our 
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theoretical  studies  led  us  to  the  concepts  of  energy  storage  in 

metastable  states  that  are  long  lived  because  of  radiation  trapping. 

We  also  proposed  to  use  an  anti-Stokes  process  to  extract  the  stored 

3 

energy  and  suggested  a  specific  process  in  atomic  mercury.  The  energy 
storage  method  and  the  anti-Stoke  extraction  process  were  subjects  of  a 
patent  application.  However,  the  application  is  now  dormant  due  to 
lack  of  Stanford  interest. 

Our  experimental  work  in  Hg^  led  almost  immediately  to  the  demon¬ 
stration  of  a  photodissociation  laser  in  Hg  which  results  from  the 

4 

dissociation  of  optically  excited  Hg^.  This  laser  also  represents  a 
new  class  of  devices  and  is  an  early  form  of  recent  photodissociation 
lasers  based  on  alkali  salts.  Of  particular  interest,  we  observed  laser 

o 

action  in  atomic  mercury  on  the  strong  5461  A  resonant  transition. 

Our  main  effort  in  these  studies  was  to  investigate  molecular  Hg^ 
as  a  potential  laser  medium. ^  Our  measurements  led  to  two  principal 
results  that  affect  the  design  of  all  high  power  gas  laser  sources.  The 
first  is  that  upper  level  absorption  in  molecular  systems  like  leads 
to  net  loss  for  low  cross-section  laser  transitions.  The  second  is  that 
three-body  collision  rates  quench  the  storage  levels  so  rapidly  that 
energy  storage  in  gases  is  not  a  useful  approach  for  high  energy  density 
storage  laser  devices. 

We  specifically  demonstrated  these  results  for  Hg^  where  we  measured 
gain,  loss,  and  collisional  deactivation  of  the  storage  level.  We  also 
stated  as  a  conclusion  of  our  work,  that  low  cross-sections  for  energy 
storage  was  not  appropriate  in  gaseous  systems.  This  point  has  now  been 
conceded  after  three  years  of  further  research  into  the  group  six  laser 
systems. 
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Our  present  understanding  of  high  power  or  energy  storage  laser 
systems  is  further  modified  by  optical  damage  problems.  High  energy 
storage  laser  systems  must  be  designed  not  to  exceed  the  damage  limit 
of  the  resonator  or  external  optical  surfaces.  For  most  laser  systems 
this  puts  an  effective  limit  on  energy  storage  density  and  thus  on  the 
minimum  useful  cross-section  required  for  efficient  energy  extraction. 
Interestingly,  Nd: Glass  is  remarkably  close  to  the  optimum  energy 
storage  laser  system  based  on  these  considerations  alone. 

I.  Two  Photon  Laser  Studies 

To  date,  population  inversion  and  subsequent  stimulated  emission 
via  the  two  photon  process  has  not  been  achieved.  The  interest  in  such 
a  system  is  the  potential  for  large  energy  storage  and  the  potential  for 
very  large  tunable  output  powers. 

We  investigated,  in  detail,  a  number  of  two  photon  laser  systems 

which  included  atomic,  molecular  and  solid  state  schemes.  Our  interests 

focused  on  atomic  systems  where  the  energy  storage  could  be  accomplished 

by  optical  pumping  and  the  two  photon  laser  action  triggered  by  a  high 

peak  power  infrared  laser  source.  Atomic  mercury,  copper,  calcium, 

strontium  and  bariam  were  theoretically  investigated.  We  also  looked  at 

two  photon  gain  in  nitrogen,  oxygen  and  thallium.  The  result  of  the 

18 

study  showed  that  the  required  energy  storage  density  was  near  10  atoms/cc 
for  at  least  1  ysec  duration.  Our  studies  further  showed  that  the  two 
photon  energy  storage  level  must  not  be  radiatively  coupled  to  lower 
levels.  The  large  energy  storage  density  requirements  led  to  a  high  pump 
laser  energy  level  of  greater  than  100  mJ.  At  the  time  of  the  study  the 
energy  was  greater  than  that  available  from  our  existing  tunable  laser 


sources. 


.1  . 


As  a  second  part  of  the  program  we  investigated  the  triggering 
requirements  to  initiate  two  photon  laser  action,  assuming  that  the 
inversion  was  available.  This  calculation,  which  is  a  major  part  of 
H.  Komine's  Ph.D.  thesis,**  showed  that  efficient  energy  extraction  from 
a  two  photon  laser  required  a  much  greater  initiation  pulse  energy  than 
was  previously  thought.  At  the  time  of  the  program  our  existing  tunable 
laser  did  not  possess  adequate  energy  for  triggering  a  two  photon  laser 
system. 

We  have  recently  reconsidered  the  two  photon  laser  source.  The 
re-evaluation  was  brought  about  by  our  Raman  studies  and  by  the  high 
pulse  energy  Nd: Glass  tunable  laser  source  now  under  construction.  The 
laser  technology  and  the  understanding  of  the  energy  storage  and  laser 
trigger  requirements  make  pursuit  of  a  two  photon  laser  a  much  more 
probable  success  now  than  four  years  ago.  The  benefits  of  high  peak 
power  tunable  output  are  still  present.  An  additional  benefit  is  the 
possibility  of  stimulated  anti-Stokes  Raman  up-conversion  for  both  fre¬ 
quency  shifting  and  energy  gain  toward  the  ultraviolet. 

J.  Frequency  Conversion  via  Four  Wave  Mixing 

We  have  recently  completed  theoretical  and  experimental  studies  of 

four  wave  mixing  as  a  frequency  conversion  process.  The  initial  concept 

was  proposed  in  1972  by  Harris  and  Byer,  but  was  not  demonstrated  until 

1977.7  The  early  experiment  showed  that  the  1.4  pm  -  4,0  pm  0P0  tunable 

source  could  be  frequency  extended  to  cover  the  4  pm  -  18  pm  spectral 

7 

range  by  four  wave  mixing  in  gas. 
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However,  the  early  experiment  was  not  in  quantitative  agreement 
with  theory  in  that  the  observed  mixing  efficiency  was  a  factor  of  30 


less  than  predicted. 

In  1978  -  1979  we  carefully  studied  the  four  wave  mixing  process 
in  gas.  During  these  measurements  great  care  was  taken  to  define 
the  spatial  and  temporal  aspects  of  the  laser  beams.  The  results  were 
the  major  part  of  S.  J.  Brosnan's  Ph.D.  thesis  which  was  completed  in 
April  1979. 8 

Our  concept  of  four  wave  mixing  and  a  method  to  implement  it  led 
to  a  patent  for  the  approach.  Further,  four  wave  mixing  was  shown  to  be 
an  integral  part  of  stimulated  Raman  scattering.  Our  increased  under¬ 
standing  of  the  interaction  of  four  wave  mixing  and  stimulated  Raman 
scattering  led  to  the  successful  demonstration  of  16  pm  generation  by 
four  wave  mixing,  assisted  stimulated  rotational  Raman  scattering  in 
para-H^  gas  using  a  CO^  laser  source. 

K.  Optical  Parametric  Oscillator  and  Optical  Parametric  Amplifier 

Studies 

The  need  for  higher  power  tunable  infrared  and  visible  sources  led 
us  to  investigate  improvements  in  the  growth  of  LiNbO^  by  increasing  the 
available  crystal  size.  Our  initial  growth  studies  showed  that  large 

9 

high  quality  LiNbO^  crystals  could  be  grown  along  the  01.4  direction. 
This  growth  direction  is  a  few  degrees  off  of  the  optimum  phasematching 
angle  for  1.06  ym  pumped  parametric  oscillation. 

We  quickly  demonstrated  that  01.4  grown  LiNbO^  could  be  used 
for  parametric  oscillation  over  an  extended  1.4  -  4,0  ym  tuning  range. 
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Our  work  then  proceeded  toward  the  development  of  a  computer  controlled 
widely  tunable,  high  energy  source. 

During  the  past  three  years  we  have  achieved  significant  progress 
toward  our  goal  of  a  fully  automatic  high  resolution,  reliable,  tunable 
source. 

We  Improved  the  LiNbO^  crystal  growth  by  continued  research  which 
led  to  a  01.3  growth  axis  and  to  crystals  20  mm  in  diameter  by  6  cm  in 

length.  These  crystals  allow  up  to  0.1  J  of  tunable  energy  to  be 

generated  at  10  pps  when  pumped  by  a  Nd:YAG  source.  We  transferred  the 
growth  technology  to  a  commercial  firm  which  now  supplies  more  than  50 
crystals  per  year  to  customers  around  the  world. 

We  developed  a  reliable  Q-switched  NdrYAG  pump  source  for  the 
0P0  measurements.  A  critical  step  in  that  work  was  the  development  of 

the  unstable  resonator  Nd:YAG  source11  which  is  now  the  basis  for  the 

standard  commercial  Nd:YAG  laser  source. 

We  developed  the  first  computer  controlled  tunable  laser  source 

with  the  capability  to  continuously  scan  over  the  widest  wavelength 

range  yet  achieved  in  a  single  device.  Recently  our  computer  software 

and  control  capability  has  been  transferred  to  commercial  products. 

We  demonstrated  continuous  tuning,  data  recording  and  data  processing  at 

10  pps  at  0.6  cm  1  linewidth  over  the  1.4  -  4.0  pm  range.  We  used  our 

12 

tunable  source  capability  in  remote  air  pollution  measurements  and  in 

13 

spectroscopic  studies.  Recently,  we  have  successfully  demonstrated 
continuous  high  resolution  tuning  at  0.08  cm  1  linewidth  over  a  20  cm  1 
range  and  have  taken  long  path  absorption  spectra  of  the  atmosphere  with 
the  source.  The  0P0  tuner  has  successfully  operated  for  longer  than  one 
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year  without  crystal  or  component  damage.  Finally,  we  have  demonstrated 
single  axial  mode  operation  of  the  0P0  and  have  fully  characterized 
crystal  damage  limitations,  optimum  OPO  design  and  optimum  pumping  con¬ 
ditions.  ^ 

Our  work  on  OPO  tuners  led  to  the  consideration  of  optical  parametric 

amplifiers  as  efficient  power  boosting  devices  with  a  wide  tuning  range. 

We  successfully  completed  a  detailed  study  of  both  a  KD*P  OPA  and  a 

LiNbO^  OPA  and  for  the  first  time  verified  the  OPA  theory  by  careful 

quantitative  measurements. There  is  no  doubt  that  OPA's  will  be  used  as 

wide  bandwidth  efficient  power  amplifiers  in  tunable  laser  systems. 

Our  work  in  this  area  was  recently  published  in  two  papers  in  the 
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Journal  of  Quantum  Electronics.  ’  The  work  formed  part  of  the  Ph.D. 
thesis  research  of  Dr.  S.  J.  Brosnan  and  Dr.  R.  Baumgartner. 

L.  Efficient  Frequency  Conversion  via  Stimulated  Raman  Scattering 

Our  work  in  four  wave  mixing  and  our  development  of  the  high  peak 
power  Nd:YAG  source  led  to  the  use  of  Stimulated  Raman  Scattering  (SRS) 
to  efficiently  convert  laser  radiation  to  new  frequencies. 

During  our  Nd:YAG  pumped  Hg^  and  OPO  studies  we  used  SRS  as  a  con¬ 
venient  means  of  efficiently  generating  new  coherent  wavelengths  at  both 
the  Stokes  and  anti-Stokes  shifted  wavelengths.  Our  work  led  to  a  better 
understanding  of  the  Raman  process^  and  to  its  general  acceptance  as  a 
useful  wavelength  shifting  method. ^  SRS  has  now  reached  commercializa- 

O 

tion  and  is  having  significant  impact  because  of  the  extended  1600  A 
to  10  Um  tuning  range  made  available  by  SRS  from  a  tunable  dye  laser 
source  pumped  via  the  unstable  resonator  Nd:YAG  source. 
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Our  early  work  in  SRS  led  to  the  consideration  of  applying  the 
technique  in  the  infrared.  A  calculation  showed  that  SRS  could  pro¬ 


vide  a  useful  method  to  efficiently  convert  10  ym  to  16  yra.  The  first 
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paper  proposing  this  concept  was  published  by  Byer  in  1976.  The 
proposal  was  met  by  general  skepticism.  However,  we  pursued  the  idea 
experimentally. 

Our  work  in  this  area  was  partially  funded  by  the  Air  Force  and 

led  to  the  concept  of  a  multiple  pass  cell  for  enhancing  the  Raman  gain. 

19 

The  paper  describing  the  cell  is  to  be  published  soon,  and  a  patent 

is  pending  under  the  Air  Force  support. 

Our  persistent  research  effort  led  to  the  first  successful  stimu- 

20 

lated  rotational  Raman  scattering  from  para-^  at  10.6  ym.  We  demon¬ 
strated  50  mj  of  16  ym  output  but  more  importantly  a  very  high  40%  photon 

21 

conversion  efficiency.  Our  results  led  to  a  careful  theoretical 

evaluation  of  the  optimum  pump  laser  source  and  to  a  proposal  for  research 

22 

in  this  area  to  AF0SR.  Stimulated  Rotational  Raman  Scattering  is  a 

conceptional  breakthrough  in  infrared  tunable  sources  in  that  it  is  not 

limited  in  peak  or  average  power  by  the  nonlinear  element  but  only  by 

the  laser  source.  It  is  readily  possible  to  design  and  demonstrate  tunable 

IR  output  in  the  important  8  ym  to  20  ym  range  at  the  Joules  per  pulse 

level  at  over  1  kHz  repetition  ratel  A  step  in  this  direction  was 

23 

recently  taken  at  Exxon  where  1.6  J  at  16  ym  was  generated  at  1  pps 
using  an  approach  based  on  our  research  and  our  multiple  pass  cell  concept. 

The  research  in  Stimulated  Rotational  Raman  Scattering  was  the 
major  part  of  W.R.  Trutna's  Ph.D.  thesis.  Dr.  Trutna  was  a  joint  student 
under  Professor  A.E.  Siegman  and  Professor  R.L.  Byer. 
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Our  work  in  stimulated  Raman  scattering  is  continuing  with  the 
further  development  of  theory  and  the  application  of  SRS  as  a  means  of 
extending  the  tuning  range  of  a  Nd: Glass  laser  source.  We  are  also 
planning  experiments  to  demonstrate  a  Raman  transfer  laser  and  efficient 
picosecond  pulse  conversion  via  the  stimulated  Raman  process. 

M.  Current  Research 

Our  research  in  laser  physics  and  laser  techniques  is  currently 
generating  new  results  in  two  areas.  We  have  recently  demonstrated 
Fourier  transform  limited  single  axial  mode  operation  of  an  unstable 
Nd:YAG  laser  oscillator  via  injection  locking.  This  work  is  part  of  the 
Ph.D.  thesis  of  Y.K.  Park  and  represents  a  technical  achievement  of  some 
importance  for  both  spectroscopic  and  military  applications  of  the  Nd:YAG 
source. 

We  are  also  pursuing  the  use  of  laser  induced  ionization  as  a  new 
and  very  selective  and  sensitive  detection  and  measurement  tool.  Research 
in  this  area  will  be  the  focus  of  our  future  research  under  this  program. 

We  plan  to  continue  to  improve  laser  sources  and  tunable  laser 
capability  in  support  of  our  research  effort.  We  also  plan  to  continue 
to  work  toward  obtaining  patents  whenever  appropriate  and  to  transfer  I 

our  laser  capability  to  the  commercial  community.  Finally  we  recognize  • 

1 

our  responsibility  to  inform  the  AFOSR  of  any  potential  benefit  to  the  [ 

Air  Force  resulting  from  our  research. 
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annual  meeting  in  San  Francisco,  October  1978. 
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APPENDIX  I 

Saturable  absorber  overlap  of  iodine  127  with  the  optically 
pumped  546-nm  mercury  lasera) 

Max  Artusy  and  A.  E.  Siegman 

Edward  L.  Gi niton  Laboratory.  Department  of  Electrical  Engineering,  Stanford  University,  Stanford. 

California  94305 

(Received  23  May  1977;  accepted  for  publication  23  June  1977) 

A  single  isolated  saturated-absorption  feature  in  l},I}  using  the  optically  pumped  346. 1-nm  Hg  laser  has 
been  observed  with  low  (~  6  MHz)  resolution.  The  ’S,— t’Pj  Hg  laser  transition  can  be  centered  exactly  on 
this  feature  by  varying  the  Nj  buffer  gas  pressure  using  the  measured  pressure  tuning  rate  of  —9.4 
MHz/Torr.  This  may  provide  an  attractive  laser  system  for  absolute  wavelength  stabilization. 

PACS  numbers:  42.65.Gv,  32.70.Jz,  42.55.Hq 


The  optically  pumped  546.  1-nm  mercury  laser1'2  is 
potentially  very  attractive  for  frequency  stabilization. 
The  absence  of  any  active  discharge  in  the  all -fused  - 
quartz  single -isotope  gain  cell  should  imply  a  nearly 
ideal  long-lived  laser  medium  especially  suitable  for 
high -stability  applications,  with  minimum  disturbance 
from  gas  cleanup,  discharge  instabilities,  and  catapho- 
resis.  Laser  oscillation  is  obtained  by  pumping  the 
continuously  flowing  or  sealed -off2  Hg-N2  laser  tube 
with  specially  constructed  low-pressure  Hg  lamps  in  a 
conventional  double -ellipse  pump  cavity.  A  simplified 
energy -level  diagram  of  the  laser  system  is  shown  in 
Fig.  1.  The  253.7-nm  Hg  resonance  line  populates  the 
3P,  level  of  mercury,  which  is  quenched  to  3P0  with  the 
aid  of  ~30  Torr  of  N2.  The  less-intense  mercury  404.7- 
nm  line  then  pumps  from  3/>0  to  the  Hg  3S,  level,  from 
which  laser  actions  occur  down  to  the  3P2  level  on  the 
mercury  green  line  at  546.1  nm.  The  nitrogen  also 
quenches  the  3P2  lower  laser  level,  without  totally 
quenching  the  other  excited  mercury  levels  involved  in 
the  laser  action. 

Intense  non -self -reversed  pump  lamps  capable  of 
operating  with  very  small  quantities  of  single  isotope 
mercury  have  been  developed  for  pumping  this  laser. 


SIMPUIVIEO  Ha  -  LEVEL  SCHEME 


FIG.  1.  Energy-level  scheme  of  the  optically  pumped  llg-N2 
laser. 


•’Research  was  supported  by  the  Air  Force  Office  of  Scientific 
Research.  Thanks  are  due  to  Nell  Holmes  who  has  carried 
out  other  aspects  of  the  Hg  laser  development  and  on  whoso 
apjjoratus  some  of  these  data  were  taken. 

Applied  Physics  Letters,  Vol.  31,  No.  5,  1  September  1977 


As  with  other  metal -vapor— rare-gas  discharges, 
cataphoresis  is  a  problem  of  primary  importance. 3  To 
overcome  this  difficulty,  we  use  the  diffusion -bypass 
lamp  design  shown  in  Fig.  2.  A  2-A  dc  discharge  in  the 
temperature-controlled  capillary  bore  of  this  lamp 
cataphoretically  pumps  mercury  front  a  cold  reservoir 
near  the  anode  to  the  cathode  at  a  rapid  rate.  The 
parallel  large -bore  discharge -free  tube  then  provides  a 
diffusion  return  path  through  which  the  mercury  can 
continuously  diffuse  back  to  the  anode.  Nickel  baffles 
suppress  the  discharge  in  the  return  path  while  allow¬ 
ing  free  diffusion  of  the  mercury.  The  vapor  pressure 
within  the  lamp  is  varied  by  heating  the  recirculating 
deionized  water  and  thus  controlling  the  capillary  bore 
temperature.  All  other  parts  of  the  lamp  are  kept 
warmer  than  the  “50  °C  water  temperature  to  prevent 
mercury  condensation  elsewhere  in  the  lamp.  It  has 
been  found  experimentally  that  the  ion -laser -type 
tungsten  dispenser  cathodes  do  not  absorb  or  getter 
the  mercury. 4  Typical  lamp  fill  is  15  mg  of  202Hg  and 
500  p  of  research-grade  argon.  Operating  lifetime  of  a 
carefully  prepared  lamp  appears  to  be  at  least  several 
hundred  hours  without  signs  of  deterioration. 

Preliminary  inverse  Lamp-dip  experiments  using 
the  mercury  laser  with  an  ,2,l2  saturable  absorption 
cell  have  been  performed.  A  iodine  cell  was  placed 
within  the  1.25-m-long  cavity  of  a  piezotuned  flowing- 
gas  mercury  laser.  The  Ilg-N2  laser  itself  has  a  pres¬ 
sure-broadening  coefficient  of  “20  MHz/Torr,4  and 
thus  presumably  operates  in  the  transition  region  be¬ 
tween  fully  homogeneously  and  inhomogeneously 


FIG.  2.  DUfusion-bypass  low-pressure  Ar-Hg  lamp  developed 
for  pumping  the  Hg-N2  laser. 
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broadened  limits.  In  agreement  with  this,  no  laser 
Lamp  dip  has  been  observed,  although  multiple  axial 
modes  will  oscillate  simultaneously.  Single  axial  mode 
operation  was  obtained  by  the  pragmatic  approach  of 
reducing  the  pumping  level  to  only  slightly  above  thresh¬ 
old,  giving  a  120 -MHz  free -tuning  range  for  the  laser 
at  fixed  gas  pressure.  Within  the  available  pressure  and 
piezoelectric  tuning  range  of  this  laser,  only  one  in¬ 
verted  Lamp  dip  was  observed,  with  a  typical  appear¬ 
ance  as  shown  in  Fig.  3.  This  single  feature  has  a 
FWHM  of  ~6  MHz,  approximately  comparable  to  I2  plus 
He-Ne  systems  previously  investigated. 5 

One  of  the  attractive  features  of  this  system  is  the 
readily  observed  pressure  shift  of  the  laser  gain  profile 
produced  by  the  N2  buffer  gas.  The  Hg  laser  can  be 
operated  over  a  N2  pressure  range  from  at  least  11  to 
40  Torr.  The  pressure  shift  of  the  laser  transition  over 
the  central  portion  of  this  range  has  been  determined 
by  observing  the  relative  position  of  the  I2  inverted 
Lamb  dip  within  the  center  of  the  laser’s  piezotuning 
range  for  various  fixed  gas  pressures.  The  results  of 
these  measurements,  as  summarized  in  Fig.  4,  show  a 
pronounced  red  shift  of  -9.4  MHz/Torr.  By  adjusting 
the  nitrogen  pressure  to  approximately  25  Torr  (which 


FIG.  4.  Frequency  shift  of  the  Hg  546.1-mm  laser  transition 
relative  to  the  mI2  saturated-absorption  peak  as  a  function  of 
N2  pressure  in  the  HG-N2  laser  tube. 


is  not  too  far  from  the  optimum  value  for  laser  opera¬ 
tion),  the  laser  gain  profile  can  be  centered  exactly 
below  the  inverted  Lamb  dip.  This  can  minimize  the 
necessity  for  third -derivative  locking  schemes  to  sup¬ 
press  slope  offset  problems  encountered  with  the  He-Ne 
and  other  lasers. 5  ^ 

j Vote  added  in  proof.  The  iodine  transition  responsible 
for  the  saturated  absorption  feature  has  tentatively 
been  identified  as  the  K(122)  line  of  the  t'"—  2  to  v'—  34 
band. 


*N.  Djeu  and  D.  Burnham,  Appl.  Phys.  Lett.  25,  350  (1974). 
2Max  Artusy,  Neil  Holmes,  andA.E.  Siegman,  Appl.  Phys. 
Lett.  28,  133  (1976). 

3K.G.  Hernqvist,  Appl.  Phys.  Lett.  16,  1  (1970). 

4Neil  C.  Holmes  (private  communication). 

HV.G.  Schweitzer,  Jr.,  E.G.  Kessler,  Jr.,  R.D.  Deslattes, 
H.P.  Layer,  and  J.R.  Whetstone,  Appl.  Opt.  12,  2927 
(1973). 
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IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS,  VOL.  QE-15,  NO.  4,  APRIL  1979 


Laser  Induced  Photoacoustic  Grating  Effects  in 

Molecular  Crystals 

JOSE  R.  SALCEDO  and  ANTHONY  E.  SIEGMAN,  fellow,  ieee 


Abstroct~\  transient  grating  of  singlet  electronic  excited  states  is 
produced  in  a  pentacene-doped  p-terphenyl  molecular  crystal  by 
optical  absorption  from  two  crossed  time-coincident  picosecond  ex¬ 
citation  pulses  at  532  nm.  The  diffraction  properties  of  this  volume 
grating  are  probed  by  a  weak,  variably  delayed,  Bragg-matched  pico¬ 
second  probe  pulse.  At  high  excitation  intensities,  a  strong  oscillatory 
behavior  in  the  time-dependent  scattering  efficiency  is  observed  super¬ 
imposed  on  the  exponential  decay  pattern  of  the  excited-state  grating. 
We  attribute  the  oscillatory  behavior  to  a  thermal  grating  which  induces 
coherent  microwave  acoustic  phonons.  These  in  turn  modulate,  at  the 
sound  frequency,  the  optical  absorption  properties  of  the  pentacene 
molecules  in  the  excited  state.  An  acoustically  induced  amplitude 
grating  effect  is  thus  obtained,  in  contrast  to  conventional  acousto¬ 
optic  phase  grating  effects. 

I.  Introduction 

E  report  the  observation  of  strong  photoacoustic  ef¬ 
fects  in  molecular  crystals,  and  particularly  the  observa¬ 
tion  of  an  acoustooptic  amplitude  grating  effect,  obtained 
while  applying  the  picosecond  transient  grating  method  to  the 
study  of  energy  migration  processes  in  molecular  crystals. 

The  details  of  our  transient-grating  method  and  of  the  en¬ 
ergy  migration  studies  have  been  given  elsewhere  [1],  [2],  but 
it  may  be  appropriate  to  briefly  review  the  transient  grating 
concept.  Two  coherent,  time-coincident  picosecond  excita¬ 
tion  pulses  are  transmitted  through  an  experimental  sample 
which  is  absorbing  at  the  excitation  wavelength,  as  in  Fig.  1. 
In  the  experiment  described  here  the  absorbers  consist  of  or¬ 
ganic  molecules  (pentacene)  imbedded  in  a  transparent  crystal¬ 
line  host  (p-terphenyl).  Optical  absorption  in  the  interference 
pattern  between  the  two  excitation  beams  produces  a  volume 
grating  pattern  of  electronic  singlet  excited  states  of  the 
pentacene  molecules.  This  excited-state  grating  then  acts  as 
an  optical  absorption  grating  or  “transient  hologram.”  The 
time-dependent  decay  of  the  excited-state  grating  is  probed 
by  Bragg  diffraction  of  a  weak,  variably  delayed  picosecond 
probe  pulse.  The  intensity  of  the  diffracted  probe  pulse  is 
monitored  as  a  function  of  probe  pulse  delay.  In  the  case  at 
hand,  the  grating  pattern  decays  through  a  combination  of 
excited-state  decay  and  excited-state  transport  (which  washes 
out  the  grating  fringes).  For  purely  diffusive  excited-state 
energy  transport,  the  decay  is  exponential,  and  the  decay  rate 
gives  the  transport  parameters  directly.  A  diffusion  coef- 
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*  DIFFRACTED 


Fig.  1.  Transient  grating  concept  and  geometry.  Probe  pulse  is  orthog¬ 
onally  polarized  relative  to  excitation  pulses,  and  signal  is  extracted 

with  cube  polarizer. 

ficient  for  the  singlet  excited  electronic  states  of  pentacene 
in  p-terphenyl  has  thus  been  obtained  [2] . 

When  we  performed  the  same  experiments  at  optical  excita¬ 
tion  intensities  above  ~500  MW/cm*,  a  strong  well  defined 
oscillatory  behavior  in  the  scattering  efficiency  versus  decay 
was  observed.  In  this  paper  we  interpret  this  oscillatory  be¬ 
havior  as  due  to  strong  photoacoustic  effects.  At  high  excita¬ 
tion  intensities,  a  significant  thermal  grating  is  superimposed 
on  the  excited-state  grating.  This  thermal  grating  produces  a 
transient  disturbance  which  can  be  described  as  a  “frozen” 
stress  pattern  plus  two  counterpropagating  coherent  micro- 
wave  acoustic  waves  [3].  The  density  modulation  associated 
with  this  disturbance  modulates,  at  the  microwave  sound  fre¬ 
quency,  predominantly  the  optical  absorption  cross  section  of 
the  pentacene  excited  electronic  singlet  states,  most  likely  by 
pressure  tuning  the  absorption  line.  This  in  turn  modulates 
the  effective  depth  of  the  excited-state  grating  and  thus  its 
scattering  efficiency,  producing  an  amplitude  grating  effect. 
To  our  knowledge,  such  photoacoustic  amplitude  grating  ef¬ 
fects  have  not  been  previously  observed.  They  provide  inter¬ 
esting  information  on  phonon  processes  and  on  pressure  tun¬ 
ing  in  the  crystals,  in  addition  to  preserving  the  excited-state 
transport  information. 

II.  Experimental  Observations 
A.  Apparatus 

The  experimental  apparatus  for  these  studies  is  outlined  in 
Fig.  2.  A  continuously  pumped,  high  repetition  rate,  Q- 
switched  and  actively  mode-locked  Nd:YAG  laser  [4]  pro¬ 
duces  bursts  of  mode-locked  pulses,  from  which  single  pulses 
are  selected  by  a  LiNbOj  Pockets  cell  in  a  Blumlein  configura¬ 
tion.  These  pubes  are  efficiently  doubled  lo  532  nm  by  a 
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Fig.  2.  Transient  grating  experimental  setup. 


temperature-tuned  90°-phase-matched  CD*A  crystal,  and  split 
into  two  strong  excitation  pulses  (==300  kW  peak  power  each) 
and  a  weak  variably  delayed  probe  pulse  (*=20  kW).  All 
pulses  are  linearly  polarized  in  the  plane  of  the  figure,  and 
are  Gaussian  in  time  and  in  space  with  a  duration  [full  width 
at  half  maximum  (FWHM)]  of  50-75  ps.  Due  to  the  par¬ 
ticular  optical  absorption  properties  of  pentacene  at  532  nm, 
excitation  and  probe  pulses  are  chosen  to  be  at  the  same  wave¬ 
length.  To  meet  the  Bragg  Condition,  the  probe  pulse  is 
brought  from  the  back  of  the  crystal  along  the  reverse  path  of 
one  of  the  excitation  pulses,  but  with  its  polarization  rotated 
by  90°.  The  first-order  Bragg-diffracted  pulse  preserves  this 
polarization  and  is  scattered  along  the  reverse  path  of  the 
other  excitation  pulse.  A  polarization  beam  splitter  then  de¬ 
flects  the  signal  into  a  photodetector  connected  to  a  lock-in 
amplifier.  The  lock-in  output  drives  they  axis  of  an  X-Y  re¬ 
corder,  while  the  x  axis  is  driven  by  the  motorized  probe  pulse 
delay  line.  For  slow  delay  scan  rate  and  high  laser  repetition 
rate  (100  Hz),  data  gathering  is  essentially  continuous  in  a  stro¬ 
boscopic  sense.  Each  experimental  plot  contains  typically 
50  000-100  000  laser  shots  and  takes  10  min  to  complete.  Ob¬ 
served  scattering  efficiencies  are  typically  several  percent,  yield¬ 
ing  high  signal  to  noise  ratio  in  the  measurements  ( S/N  «  1 00  is 
easily  achieved).  To  further  enhance  the  signal  discrimination, 
every  other  pair  of  excitation  pulses  is  chopped  and  phase 
sensitive  detection  is  used.  A  microscope  objective  can  also 
be  used  to  image  the  interaction  volume  onto  a  vidicon  tube 
with  large  magnification,  allowing  the  display  of  the  grating 
fringes  on  a  television  monitor.  This  provides  a  very  accurate 
tool  to  align  all  three  beams  in  the  crystal. 

H.  Results 

Results  showing  diffusive  energy  transport  at  low  excitation 
intensity  have  been  previously  reported  and  studied  |2|.  Here 
wc  concentrate  on  the  high  excitation  intensity  results,  which 
show  the  oscillatory  character  of  the  time-dependent  scatter¬ 
ing  efficiency.  Fig.  3  shows  this  oscillatory  dependence  as  ob¬ 


oe  lay 


Fig.  3.  High  excitation  intensity  transient  gxating  results  showing  the 
dependence  of  the  oscillation  period  T  on  the  grating  fringe  spacing 
A,  for  four  experimental  results. 

served  in  four  experiments  from  a  series  performed  under  iden¬ 
tical  conditions  except  for  different  fringe  spacings  A,  which 
are  adjusted  by  changing  the  angle  between  the  two  excitation 
beams.  In  all  these  results  the  grating  k  vector  is  parallel  to 
the  crystal  b  axis.  The  time  period  T  of  the  oscillatory  modu¬ 
lation  changes  with  fringe  spacing  A  with  the  linear  depen¬ 
dence  shown  in  Fig.  4.  A  straight-line  fit  yields  a  slope  of 
2.63  X  10s  cm/s,  which  we  take  to  be  the  velocity  of  sound. 
The  velocity  of  sound  in  the  same  crystal  along  the  same  di¬ 
rection  was  independently  measured  to  be  vs  =  2.65  ±  0.05 
cm/s,  confirming  the  acoustic  origin  of  the  oscillations. 

Further  examination  of  the  shape  of  the  oscillatory  decays 
in  Fig.  3  leads  us  to  conclude  that  the  observed  effects  must 
result  from  thermal  excitation  of  an  acoustooptic  amplitude 
grating.  This  grating  must  be  produced,  we  believe, by  acoustic 
pressure  tuning  of  the  excited  molecular  state  absorption  in 
the  sample.  The  theoretical  analysis  supporting  these  con¬ 
clusions  is  outlined  in  the  following  section. 

HI.  Theoretical  Formulation 

A.  Time  Dependence  of  the  Scattering  Efficiency 

The  geometrical  and  parametric  dependence  of  the  scattering 
efficiency  of  volume  holograms  is  well  known  |5|.  |(»|.  Here 
we  arc  mainly  interested  in  the  exponentially  decaying  and 
oscillating  time  dependence  of  the  observed  scattemig  from 
the  volume  grating  patterns  produced  by  thermal  plus  excited- 
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Tin.) 

Fig.  4.  Lineai  dependence  of  A  versus  T. 


state  gratings  in  our  experiments.  In  these  experiments 
scattered  probe  signals  were  observed  primarily  at  the  First- 
order  Bragg-matching  angle.  However,  the  volume  grating 
pattern  also  contains  higher  spatial  harmonics,  particularly 
at  the  second-order  Bragg  angle,  which  should  also  be  observ¬ 
able.  In  general  the  diffracted  probe  signal  intensity  /„(/)  at 
the  nth  order  Bragg  matching  angle  is  given  by 


0l/,(O*X„(OlJ  dt' 


(1) 


where  the  *  denotes  convolution,  and  lp(t'  -  t )  and /,,(/')  are 
the  probe  and  excitation  pulse  intensities,  with  the  probe 
pulse  being  delayed  in  time  by  an  amount  /.  Out  of  the  total 
complex  susceptibility  x(x,  t)  in  the  sample  the  relevant  part 
is  the  spatial  Fourier  component  of  order  n,  i.e.,  the  compo¬ 
nent  Xn(OexP  Unkx),  which  corresponds  to  spatial  vector 
nk  =  n2rr/A,  where  A  is  the  grating  spacing.  As  written  here, 
X(x,f)  is  the  Green’s  function  or  the  unit  impulse  response 
produced  by  an  excitation  pulse  that  is  a  unit  impulse  in  time. 
In  this  paper  we  will,  in  fact,  treat  the  excitation  pulse  Ie(t) 
as  a  unit  impulse,  since  the  ~70  ps  excitation  (and  probe) 
pulses  are  short  compared  to  the  ~  1  ns  oscillation  and  decay 
times  in  these  experiments.  The  scattering  susceptibility  x«(0 
is  in  general  complex,  i.e.,  x«  =  Xn  +/X«.  with  the  real  part 
being  responsible  for  phase  grating  or  index  of  refraction 
grating  effects,  and  the  imaginary  part  being  responsible  for 
amplitude  or  absorption  grating  effects. 

The  physical  process  that  produces  the  coherent  acoustic 
waves  of  interest  here  is  a  thermally  induced  acoustic  grating 
which  is  superimposed  on  the  excited-state  grating  [2].  The 
thermal  grating  and  the  associated  acoustic  waves  can  be 
analyzed  starting  from  the  linearized  hydrodynamic  equations 
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Here  A T  and  Ap  arc  the  local  temperature  and  density  changes, 
p o  is  the  equilibrium  density,  r?  is  the  viscosity, 0  is  the  expan¬ 
sion  coefficient,  X  is  the  thermal  conductivity,  a  is  the  optical 
absorption  coefficient,  and  l.c  is  the  local  optical  electric  field 
in  the  interaction  volume.  Fleet roslrict ion  can  be  neglected, 
and  the  transient  solution  for  the  density  variation  Ap  can 
then  be  readily  obtained  |8)  foi  an  excitation  pulse  very  shoit 


compared  to  the  acoustic  times  of  interest: 

Ap(jr,  /)=  Apmax  [e  ,/r°c  cos  cuf  -  e  1 /r,h )  cos(2*a:/A). 

(3) 

Here  rac  is  the  acoustic  wave  attenuation  time  or  sound  absorp¬ 
tion  time,  Tth  is  the  thermal  conduction  time,  and  u>  =  2rr/xA 
is  the  acoustic  frequency  where  rA  =  A /us,  with  ,,  being  the 
acoustic  wave  velocity.  For  the  acoustic  frequencies  generated 
in  these  experiments  (<o/2rr  ~  300-800  MHz),  the  acoustic 
decay  time  is  Tac~  50-100  ns,  and  the  thermal  conduction 
decay  time  can  be  ignored. 

A  more  complete  analysis  would  take  into  account  the 
Gaussian  transverse  spatial  dependence  of  the  excitation 
beams,  as  well  as  the  finite  duration  of  the  excitation  pulse 
le(t).  In  particular  there  will  be  an  additional  effective  decay 
of  the  acoustic  signals  as  the  finite  trains  of  left-  and  right¬ 
going  acoustic  cycles  move  across  each  other  and  eventually 
cease  to  overlap  spatially  with  each  other.  Since  our  gratings 
typically  had  at  least  20  fringes  across  the  half  width  of  the 
beam  spot  size,  however,  this  had  little  effect  during  the  5  or 
6  acoustic  cycles  typically  followed  in  our  experiments. 

Although  we  will  speak  loosely  of  pressure  tuning  effects, 
we  attribute  the  observed  effects  in  our  experiments  to  in¬ 
duced  changes  in  the  local  optical  properties  of  the  medium 
caused  primarily  by  local  density  changes  rather  than  either 
local  temperature,  local  stress,  or  pressure  effects.  Fig.  5  is 
a  plot  of  the  excited  state  grating  pattern  Nx  (x,  t  -  0),  the 
ground  state  molecular  density  N0(x,  t  =  0),  and  the  instanta¬ 
neous  local  density  change  A p(x,t),  with  the  latter  quantity 
being  plotted  at  successive  instants  of  time  through  one  com¬ 
plete  acoustic  cycle,  neglecting  damping  or  decay  effects.  The 
acoustic  density  grating  pattern  is  in  phase  spatially  with  the 
excited-state  grating  pattern,  as  well  as  with  the  sinusoidal  part 
of  the  ground  state  molecular  density.  Because  of  the  instan¬ 
taneous  nature  of  the  thermal  excitation  in  bulk,  both  the  den¬ 
sity  and  the  pressure  must  be  described  as  having  a  static  or 
“frozen”  grating  component  as  well  as  two  counterpropagating 
traveling-wave  components.  (The  “frozen”  components  will  of 
course  decay  away  with  3  much  longer  thermal  relaxation  time.) 
Note  also  that  the  density  variations  with  time  at  any  one  point 
in  space  arc  “one-sided,”  i.e„  the  local  density  swings  up  and 
back,  or  down  and  back,  but  not  both  up  and  down  at  any  one 
location.  The  molecular  densities  N0(x)  and  A, (a)  decay 
slowly  due  to  cxcited-state  relaxation  and  diffusion. 

B.  Phase  Versus  Amplitude  Grating  Effects 

Conventional  acoustooplic  interactions  in  solids  neaily  al¬ 
ways  occur  through  a  phase  grating  mechanism,  caused  by  a 
change  in  local  index  of  refraction  associated  with  the  acoustic 
disturbance.  Examination  of  the  theoretical  expressions  and 
of  Fig.  5  shows,  however,  that  this  cannot  be  the  primary  ex¬ 
planation  for  our  observations,  and  that  we  are  clearly  observ¬ 
ing  an  acoustooptic  amplitude  grating  effect. 

In  our  results  the  instantaneous  scattering  efficiency  is  seen 
to  be  reduced  nearly  to  zero  at  each  odd  half  acoustic  cycle 
after  the  excitation  pulse.  This  is  also  the  instant  when  the 
density  grating  reaches  its  maximum  amplitude  in  spatial 
phase  with  the  cxcited-state  grating.  The  implication  is  that 
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Fig.  5.  Density  wave  dynamical  behavior.  Fringe  peaks  are  atrc/A  =  0, 
1,  2,  •  ■  - .  Excited  and  ground  state  population  spatial  distribution  at 
t=  0. 

the  density-induced  scattering  at  those  instants  must  add  in 
opposite  phase  (not  in  quadrature)  to  the  excited-state  in¬ 
duced  grating.  We  believe  the  excited-state  grating  must 
surely  be  an  amplitude  grating,  and  hence  the  density-induced 
grating  must  be  an  amplitude  rather  than  a  phase  grating  also. 

Equation  (3)  shows  that  the  density  grating  consists  of  two 
traveling-wave  components  of  relative  amplitude  each,  and 
a  fixed  or  constant  component  of  relative  amplitude  - 1 .  The 
two  traveling-wave  components  decay  with  the  acoustic  decay 
time  Tac,  which  lasts  25-50  acoustic  cycles,  while  the  fixed 
component  decays  with  the  thermal  decay  time  rth ,  which  is 
much  longer.  Each  of  these  components  will  induce  a  scat¬ 
tered  signal  field  component  at  the  Bragg  angle,  with  cor¬ 
responding  amplitudes  and  phases.  Note  that  the  three  den¬ 
sity  grating  components,  and  also  the  three  vectorial  scattered 
wave  components,  sum  to  exactly  0  at  t  =  0  since  the  density 
modulation  starts  from  an  initial  value  that  is  identically  zero. 

Fig.  6  shows  vectorially  in  the  complex  plane  the  scattered 
vector  field  amplitudes  contributed  by  each  of  these  three  den¬ 
sity  wave  components,  plus  the  larger  component  contributed 
by  the  excited-state  grating,  at  three  instants  of  time  corre¬ 
sponding  to  0,  j,  and  I  acoustic  cycle.  The  sketches  are 
drawn  with  vertical  arrows  representing  the  excited  state  grat¬ 
ing,  assumed  to  be  a  pure  amplitude  grating.  The  heavy  arrow 
in  each  case  is  the  vector  sum,  or  the  instantaneous  vector 
amplitude  of  the  total  Bragg-scattered  signal  field.  The  length 
of  this  vector  squared  would  give  the  total  Bragg-scatlcrcd  in¬ 
tensity  at  that  instant. 

The  upper  sketch  assumes  that  the  density  waves  induce  a 
phase  grating  effect,  i.c„  the  three  acoustically  induced  com- 


AMPLITUDE  PLUS  PHASE  GRATING 


AMPLITUDE  PLUS  AMPLITUDE  GRATING 


Fig.  6.  Graphical  (phasor)  description  of  acoustically  induced  phase  or 
amplitude  grating  added  to  excited  state  amplitude  grating.  Resultant 
time  dependences  are  radically  different.  Only  acoustically  induced 
amplitude  grating  effects  can  explain  observed  time  dependence. 

ponents  are  90°  out  of  phase  with  the  excited-state  grating 
component;  while  the  lower  sketch  assumes  the  density  wave 
components  are  in  phase  with  the  excited-state  component. 
The  two  smaller  density  wave  components  rotate  in  opposite 
directions  in  the  complex  plane  at  the  acoustic  frequency, 
since  they  represent  oscillating  or  traveling-wave  terms.  Note 
that  the  observed  amplitude  variation  of  the  total  scattered 
signal  with  time  can  only  be  explained  by  the  lower  sketch. 
An  acoustic  phase  grating  would,  in  fact,  cause  the  total 
scattered  intensity  to  have  an  amplitude  variation  of  opposite 
sign. 

C.  Detailed  Analysis  of  the  Scattering  Efficiency 

Leaving  aside  inessential  constants  and  geometrical  factors, 
the  scattering  susceptibility  xB(0  can  be  written  as 

Xn(0  ~i  jZ^/.pC*.  0  l)|  (4) 

where  yV/  p  and  aip  are  the  density-modulated  molecular  state 
populations  and  cross  sections,  respectively,  and  where  the 
subscript  n  indicates  that  only  the  exp(jn2irx/ A)  spatial  com¬ 
ponent  is  to  be  considered.  The  sum  in  the  present  case  ex¬ 
tends  over  two  levels,  the  ground  state  i  =  0  and  the  first  ex¬ 
cited  singlet  level  /  =  1 ,  The  population  and  cross  sections  may 
be  written  as 

Mi.p(x,t)  =  N,(x,t)[\  +  Ap(x,f)/p0] 

°i,p(xJ)  =  Oi  +  (dOjldp)  Ap(x,t)  (5) 

where  /Vj(.v,  f)  and  o,  without  arguments  give  the  unperturbed 
values.  The  values  both  of  o(-  and  of  their  derivatives  may  in 
general  be  complex,  with  the  real  part  corresponding  to  ab¬ 
sorption  and  the  imaginary  part  corresponding  to  phase  shift 
or  index  of  refraction  effects. 

The  general  expression  for  Ap(x,r)  has  been  given  in  (3). 
The  populations  of  the  ground  and  excited  levels  produced 
by  the  excitation  beam  delta  functions  arc  [2) 
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{  Nrx  e-,lT  [1  -  e-k'Dt  cos(2»rx/A)] 

N0(x,t)  =  N0-  Ni(x,t)  (6) 

where  N0  without  arguments  is  the  total  density  of  absorbing 
molecules  in  the  crystal,  A,/na*  is  the  initial  optically  induced 
excited-state  density  at  the  grating  peaks,  r  is  the  excited 
state  fluorescent  lifetime,  and  (1/t  +  k2D)  =  K/2  is  the  decay 
rate  for  the  excited-state  fringe  pattern  due  to  excited-state 
relaxation  plus  excited-state  diffusion  with  diffusion  coef¬ 
ficient  D. 

We  now  define  the  normalized  maximum  density  variation 
n  =  Apmax/p0  «  1  (7) 

and  the  oscillatory  modulation  component 

M(t)  =  exp  (-  tlTac)  cos  Cjjt  -  1 .  (8) 

We  also  define  the  cross  section  sensitivity  parameters,  or 
normalized  density  derivatives 

50  =(9oo/dp)Po/o0 

51  =(aoi/9p)p0/oi  (9) 

which  could  be  complex  or  even  imaginary  if  the  density 
modulation  produces  an  index. of  refraction  grating.  We 
then  obtain  the  n  -  1  Bragg  angle  component  of  the  scatter¬ 
ing  susceptibility  from  (4),  (8),  and  (9)  as 

-/X.(0  =  (^m,x/2)(a,  -Oo)e<'l”k'D" 

+  Tjo,  (S,  +  1)  (Nrxl2)  e',/r  M(t) 

+  r?a0(50  +  1)  (N0  -  (7V,max/ 2)  e''/T)M(/).  (10) 

This  expression  contains  three  terms.  The  first  is  the  basic 
excited-state  transient-grating  signal  [2]  unperturbed  by 
acoustic  effects.  It  gives  the  pure  exponential  decay  observed 
at  low  intensities,  as  determined  by  the  excited-state  lifetime 
t  and  the  diffusion  coefficient  D.  The  second  and  third  terms 
describe  how  the  first  excited  state  and  ground-state  diffrac¬ 
tion  effects  are  modulated  by  the  thermally  induced  densit / 
variations  or  coherent  acoustic  waves.  In  the  (St  +  1)  factors 
in  the  second  and  third  terms,  the  factor  of  1  represents  simple 
density  modulation  of  the  local  molecular  density,  while  the 
5/  factors  represent  density  modulation  of  the  molecular  cross 
sections  a,  ,  as  in  (5).  It  will  turn  out  that  the  magnitudes  of 
the  Sj  seem  to  be  substantially  larger  than  unity  to  fit  the  ob¬ 
served  results. 

The  resulting  diffracted  signal  intensity  is  then  given  by 

/i(0*  f  t)\e-(KI2)t' 

BO 

+  A,e~' /T  M(t') 

+  /10(2A0//Vfnax  ^  c*''/T)A/(f')  I2  dt'  (11) 

where  we  have  normalized  the  initial  scattering  efficiency  to 
unity.  The  oscillatory  terms  have  magnitudes 

.  »7  (S„  +  I )  a„ 


+  l)o, 

A  i  =  - . 


Note  that  squaring  this  expression  not  only  squares  each  of  the 
above  terms,  but  most  importantly  leads  to  beating  effects 
between  the  basic  excited-state  grating  and  the  two  acoustic 
terms.  Because  the  probe  pulsewidth  of  ~70  ps  is  short  com¬ 
pared  to  the  time  constants  involved  in  Xr(0>  we  can  simplify 
(1 1)  as 

/,(')=  + 

+  /t0(2jV0/rVimax  -  e~^T)Af(/)|2 
=  *  +  AiM(t)c't/r 

+  Aoi2N0/Nlnax  -  e',/T)Af(0<?'(*/2)f-  (13) 

Here  we  have  also  ignored  products  of  order  172  in  the  second 
line,  since  tj2  «  1  • 

Fig.  7  shows  (13)  plotted  for  typical  values  of  the  relevant 
time  constants  as  we  slowly  “turn  on”  the  excited-state 
density-modulation  term  At  and  the  ground-state  density 
modulation  term  Ao  independently.  The  experimental  param¬ 
eters  employed  are 

A'rax/Aro  =  0.3 

t  =  9.5  ns 
rac  =  1 00  ns 
A  =  5  pm 

vs  =  2.65  X  105  cm/s 

D  -  0.5  cm2/s.  (14) 

It  clearly  seems  that  the  shape  of  our  experimental  results  can 
be  matched  only  by  assuming  that  density  modulation  of  the 
excited  state  is  the  dominant  mechanism  causing  the  acoustic 
effects.  This  is  supported  by  the  data  shown  in  Fig.  8,  ob¬ 
tained  from  the  third  plot  from  the  top  in  Fig.  3.  The  plot 
shows  both  the  exponential  decay  of  the  peak  amplitudes, 
which  gives  the  diffusion  coefficient  D,  and  also  the  decay  of 
the  dcpth-of-modulation  amplitudes  i.e.,  the  difference  be¬ 
tween  the  fitting  exponential  and  the  observed  amplitude  at 
the  minima.  The  decay  rates  are  essentially  the  same,  which 
says  that  the  oscillatory  term  decays  at  the  same  rate  as  the 
grating  component  of  Nt(x),  not  N0(x).  Only  if  density 
modulation  of  the  excited  state  is  the  dominant  mechanism 
can  this  depth  of  modulation  decay  be  explained.  For  ground 
state  modulation  the  decay  would  have  to  be  significantly 
slower. 

D.  Comparison  to  Experimental  Results 

Fig.  9  repeats  the  experimental  result  of  Fig.  3,  second 
curve  from  the  top  (A  =  5.33  pm),  along  with  a  semilog  plot 
of  the  peak  amplitudes  in  this  curve,  showing  their  almost 
exact  exponential  decay.  At  the  peak  times,  the  instantaneous 
density  change  Ap(x,t)  is  zero,  and  hence  the  obscivcd  expo¬ 
nential  decay  rate  for  the  peaks  should  be  (2 It  +  2 k2 1)).  The 
observed  results  are  in  fact  in  agreement  with  low-intensity, 
nonoscillating  measurements  which  give  t  =  9.5  ns  and  D  -  0.5 
cni2/s  for  this  particular  sample.  Also  shown  in  Fig.  is  a 
theoretical  fit  using  (13)  for  these  experimental  conditions, 
with  Ao  =  0.02  and  /I,  =  +0.37.  l-vidcntly  it  is  density  mod¬ 
ulation  of  the  excited  rather  than  the  ground  level  tint  is 
dominant. 
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Fig.  7.  Excited  state  (upper)  modulation  and  ground  state  (lower) 
modulation,  for  several  modulation  "depths”  A i  and /40.  Although 
behavior  may  look  similar,  the  oscillation  decay  is  radically  different, 
being  much  slower  for  ground  state  modulation. 

From  all  our  other  measurements,  we  can  estimate  that 
ffi/o0  *5,  and  from  the  heat  deposition  and  thermoacoustic 
properties  of  the  material,  we  can  estimate  that  in  the  high  in¬ 
tensity  experiments,  the  normalized  maximum  density  change 
corresponds  to  r?~10'3.  From  (12)  we  then  obtain  the 
rough  values 

(S0  +  1)  — 100 

(St  +l)~+500 


Fig.  8.  »-Pealt  amplitudes  for  third  plot  in  tig.  3  (A  =  4.7  um).  show¬ 
ing  excellent  exponential  decay.  Decay  time  constant  gives  D  -  0.S 
cm2/s  for  this  particular  sample,  for  diffusion  along  the  h  axis. 
X-Depth  of  modulation  points  obtained  from  minima  of  same 
plot.  Observed  decay  can  only  be  explained  if  an  acoustic  modula¬ 
tion  of  excited  state  optical  properties  is  assumed.  Ground  state 
modulation  would  produce  a  much  slower  decay. 


WtA 

Fig.  9.  Theoretical  fit  to  second  experimental  plot  in  Fig.  3  (A  *  5.3 
rim),  for  Aq  ■  -0.02  and  At  =+0.37. 
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(15) 


at  the  peak  of  the  acoustic  wave  (where  A p/p0  <  0).  There  is 
independent  support  for  the  assumption  that  8oo/dp<0,  if 
we  consider  the  ground  state  absorption  spectrum  of  penta- 
cene  in  the  vicinity  of  532  nni,  and  a  static  pressure  tuning 
study  of  pentaccne  lines  [10].  No  data  on  the  pressure  tuning 
of  the  first  singlet  excited  state  of  pentaccne  seem  to  be  avail¬ 
able,  but  our  results  imply  that  the  net  effect  at  532  nm  must 
be  of  opposite  sign,  and  several  times  larger.  The  agreement 
shown  in  Fig.  9  is  good,  especially  if  we  look  at  the  simplicity 
of  the  theoretical  model,  in  which  we  did  not  include  highly 
nonlinear  effects  potentially  present  in  the  interaction  volume, 
such  as  the  following. 

1)  In  the  high  compressibility  .egime  produced  by  the  laser 
induced  acoustic  waves  (excess  pressure  -  50  atm),  the  sound 
velocity  can  no  longer  be  considered  constant  [1 1|,  and  the 
acoustic  waves  can  no  longer  lie  rigorously  treated  by  a  linear- 

teed  model. 


2)  A  linear  dependence  of  cross  section  on  density  has  been 
assumed,  when  in  more  rigor  we  should  use  the  exact  nonlin¬ 
ear  absorption  lineshape  around  532  nm. 

However,  we  believe  including  these  nonlinear  effects  would 
only  change  the  fit  in  Fig.  9  by  a  negligible  amount,  and  in  the 
right  direction  for  an  even  more  exact  fit. 

IV.  Conclusion 

The  present  paper  reports  the  observation  of  amplitude  grat¬ 
ing  effects  produced  by  coherent  acoustic  waves  induced  by 
optical  absorption  of  picosecond  pulses.  These  pulses  form 
overlapping  transient  excited-state  and  thermal  gratings  in  the 
bulk  of  a  molecular  crystal.  A  simple  model  introduced  to 
treat  the  coherent  acoustooptic  interaction  fits  the  experi¬ 
mental  results  with  quite  good  agreement.  The  results  yield 
information  about  a  variety  of  processes,  including: 

1 )  The  nature  ami  rate  of  energy  migration  for  the  molec¬ 
ular  excited  states,  obtained  by  studying  the  decay  of  '.he 
oscillatory  signal  peaks. 

2)  Information  regarding  high-frequency  acoustic  plumon 
processes  in  the  crystal,  including  acoustic  velocity  and  atten¬ 
uation.  By  changing  the  grating  orientation  in  the  crystal, 


-  37  - 


256 


IEEE  JOURNAL  Ol  QUANTUM  ELECTRONICS.  VOL.  QE-15,  NO.  4,  AI’RIL  1979 


anisotropic  and  nonlinear  acoustic  effects  might  also  be 
directly  studied. 

3)  Information  regarding  the  pressure  or  density  modulation 
of  the  ground-state  and  excited-state  absorption  in  the  crystal. 

The  amount  of  information  simultaneously  obtained  from 
measurements  such  as  these,  concerning  both  molecular  and 
crystalline  properties,  is  sizable.  We  hope  this  process  can  be 
naturally  extended  to  other  materials  as  well,  providing  a 
unified  understanding  of  their  properties. 

Acknowledgment 

Earlier  work  on  the  transient  grating  experiments  was  sup¬ 
ported  by  the  Joint  Services  Electronics  Program  at  Stanford 
University  and  by  the  National  Science  Foundation.  This 
work  grew  out  of  studies  of  energy  transport  in  pentacene 
carried  out  in  collaboration  with  Prof.  M.  Fayer  and  D.  Dlott 
of  the  Stanford  Chemistry  Department,  and  we  appreciate 
the  use  of  molecular  crystal  samples  provided  by  them. 

References 

|1]  D.  W.  Philiion,  D.  J.  Kuizenga,  and  A.  E.  Sicgnian,  “Subnanosec¬ 
ond  relaxation  time  measurements  using  a  transient  induced 
grating  method,”  Appl.  Phy>.  Lett.,  vol.  27,  pp.  85-87,  1975. 


( 2 ]  J.  R.  Salcedo,  A.  E.  Sicgman,  D.  D.  Dlott,  andM.  D.  I'aycr,  ''Dy¬ 
namics  of  energy  transport  in  molecular  crystals:  The  picosecond 
transient  grating  method," Phys.  Rev.  Lett.,  vol.  41,  pp.  131-133, 
July  10,  1978;  also  in  Picosecond  Phenomena ,C .  V.  Shank,  E.  P. 
Ippen,  and  S.  L.  Shapiro,  Eds.  Berlin:  Springer  Series  in  Chcm. 
Phys.,  vol.  4,  1978,  pp.  240-243. 

(3(  D.  C.  Auth,  "New  high-power  source  of  Coherent  microwave 
phonons, "/tpp/.  Phys.  Lett. ,vol.  76, pp.  52 1-5 23, June  15, 1970. 

(4)  D.  W.  Philiion,  “Evolution  of  mode-locked  pulses  and  their  ap¬ 
plication  to  time-resolved  spectroscopy,”  Ph.D.  dissertation, 
Stanford  Univ.,  Stanford,  CA,  Nov.  1974. 

[5]  H.  Kogelnik,  “Coupled  wave  theory  for  thick  hologram  grat¬ 
ings,”^//  Syst.  Tech.  J .,  vol.  48,  pp.  2909-2947, 1969. 

(6}  A.  E.  Siegman,  “Bragg  diffraction  of  a  Caussian  beam  by  a 
crossed-Gaussian  volume  grating,”  J.  Opt.  Soc.  After.,  vol.  67, 
pp.  545-550,  1977. 

(7)  D.  R.  Dean,  “Optically  induced  diffraction  gratings  in  liquids  and 
solids,”  Office  of  Naval  Research,  Washington,  DC,  Tech.  Rep. 
TR-73-02. 

(8)  M.  E.  Mack,  “Stimulated  thermal  Rayleigh  scattering  with  pico¬ 
second  pulses, "Ann.  A.  Y.Acad.  Sc/.,  vol.  168,pp.419-436, 1970. 

(9)  R.  Herman  and  M.  Gray,  “Theoretical  prediction  of  the  stimu¬ 
lated  thermal  Rayleigh  scattering  in  liquids,”  Phys  Rev.  Lett., 
vol.  19,  pp.  824-827,  1967. 

(10)  J.  M.  Donnini,  “Effect  d’une  compression  sur  les  spectres  de 
l’anthraccne  et  du  pentacene  en  solution  dans  le  p-terphenyle 
a  1.7  K ,”/.  Own.  Phys.,  II,  vol.  11-12,  pp.  1542-1545,  1974. 

(11(  H.  Eichlcr  and  H.  Stahi,  “Time  and  frequency  behavior  of  sound 
waves  thermally  induced  by  modulated  laser  pulses,"  J.  Appl. 
Phys.,  vol.  44,  pp.  3429-3435,  Aug.  1973. 


-  38  - 


APPENDIX  III 


932' 


IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS.  VOL.  QE-IS.  NO.  9.  SEPTEMBER  1979 


Diffraction-Biased  Unstable  Ring  Resonators  with 
Possible  Applications  in  Laser  Gyroscopes 

SHINAN-CHUR  SHENG 


Abstract-The  performance  of  laser  gyroscopes  is  seriously  limited  at 
present  by  the  lock-in  phenomenon  that  occurs  between  the  counter- 
propagating  waves  at  low  rotation  rates  or  low  values  of  the  beat  fre¬ 
quency  between  the  two  waves.  To  minimize  or  eliminate  this  lock-in 
phenomenon,  different  techniques  such  as  mechanical  dither  and 
Faraday  cell  nonreciprocal-phase  elements  have  been  proposed  and 
Studied.  Here  a  new  idea  is  presented,  which  shows  how  nonreciprocal 
propagation  and  a  zero-rotation -rate  frequency  offset  can  be  obtained 
by  a  “diffraction-bias”  technique  employing  differential  spatial  filtering 
in  a  negative-branch  unstable  ring  resonator. 

I.  Introduction 

E  PROPOSE  here  a  basic  new  technique  for  obtaining  a 
zero-rotation-rate  frequency  offset  between  the  counter- 
propagating  waves  in  a  ring  laser  by  using  a  “diffraction-bias” 
technique  in  a  negative-branch  unstable  resonator.  Whether 
this  will  represent  a  useful  technique  remains  to  be  seen. 
However,  it  does  appear  to  provide  a  fundamentally  new 
attack  on  the  important  problem  of  providing  a  frequency 
offset,  and  eliminating  the  lock-in  problem  for  ring  laser 
gyroscopes  [1]  at  very  low  rotation  rates. 

The  general  structure  of  a  negative  branch  unstable  ring 
resonator  with  a  separate  focus  for  the  mode  propagating  in 
each  direction  is  shown  in  Fig.  1.  To  explain  the  basic  princi- 

Manuscript  received  March  S.  1979;  revised  M.nch  15.  1979  1  his 
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The  author  is  with  I  dwatd  l..  (bn/ ton  Laboratory.  Stanford  Univer¬ 
sity,  Stanford.  ('A  94305 
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Fig.  I.  A  typical  negative-branch  unstable  ring  resonator. 

pie  of  nonrcciprocal-diffraction  filtering,  ab  unstable  resona¬ 
tor  with  cylindrical  symmetry  along  the  optical  path  is  as¬ 
sumed.  Discussion  of  the  problem  of  astigmatism  is  deferred 
to  later  in  this  paper. 

In  Fig.  1  the  clockwise  (CW)  propagation  wave  has  a  focus 
at  f'cw,  while  the  counter-clockwise  (CCW)  wave  has  one  at 
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Fig.  2.  Schematic  representation  of  transverse  dimensions  of  CW  and 
CCW  waves  at  the  CW  focal  plane  Fqv/- 


Fig.  3.  illustration  of  differential-spatial  filtering,  forcing  oscillation 
in  l  =  1  for  the  CCW  wave,  while  allowing  1  =  0  oscillations  for  the  CW 
wave. 


Fccw-  1°  general  Few  and  Fee w  do  not  have  the  same  loca¬ 
tion  in  a  negative-branch  ring  unstable  cavity.  If  proper  de¬ 
sign  is  used  (e.g.,  equivalent  Fresnel  number  jVeq  4*0.5,  1.5, 
2.5  ■  •  • ,  where  the  higher  order  transverse  modes  are  well 
separated,  plus  appropriate  intracavity  apertures,  etc.),  the  CW 
and  CCW  waves  will  both  oscillate  in  the  lowest  transverse 
mode,  with  azimuthal  symmetry  parameter  l  =  0, 

At  the  focal  plane  containing  the  focus  Few  the  transverse 
dimensions  of  the  CW  and  CCW  waves  are  different,  as  shown 
in  Fig.  2.  At  the  other  focal  plane  containing  FCcw.  the  situa¬ 
tion  is  just  reversed,  with  the  CCW  wave  small  at  the  center 
and  the  CW  wave  larger. 

Suppose  we  insert  two  fine  pins  at  the  focal  plane  Fcw  as 
shown  in  Fig.  3  in  such  a  manner  that  they  are  not  deep 
enough  to  disturb  the  CW  wave  propagation,  but  they  do 
penetrate  into  the  region  that  the  CCW  wave  occupies.  This 
“differential-spatial  filtering”  will  cause  the  CCW  wave  to  suf¬ 
fer  very  high  diffraction  loss  in  the  /  =  0  mode,  and  hence  the 
allowed  oscillations  in  that  direction  can  only  be  the  modes 
with  /  =  1 , 2,  3,  •  •  •. 

Usually,  higher  /  azimuthal  modes  have  larger  transverse  di¬ 
mensions.  So,  an  aperture  at  the  focal  plane  Fcw  can  restrain 
the  CCW  wave  from  oscillating  in  l>  2  azimuthal  modes.  Sim¬ 
ilarly  a  proper  aperture  at  the  focal  plane  FC(  W  can  stop  any 
1 2*  1  oscillation  of  the  CW  wave. 

If  instead  of  two  pins  we  put  in  4,  6,  8  •  •  •  pins  at  equal  an¬ 
gles'  wc  may  be  able  to  force  the  lowest  azimuthal  mode  of 
the  CCW  wave  to  run  in  /  =  2, 3, 4  •  •  -,  respectively. 

Similarly,  a  pin  diffraction  filter  can  be  placed  at  focal  plane 
Fee w  to  obtain  forced  higher  azimuthal  mode  oscillation  of 
the  CW  wave.  If  CW  and  CCW  waves  arc  operating  in  different 
/  modes,  then  the  roundtrip  phase  shills  will  be  nonreciprocal 
in  the  two  directions,  and  a  frequency  offset  of  “dil fraction- 

'"Diffraction  bias"  is  not  limited  to  azimuthal  filtering  only.  In  fact, 
any  selective  filtering  mechanism  which  differentially  enhances  cer¬ 
tain  transverse  modes  (radial  or  azuniiMi.il,  .»  or  rt  m  one  dues  lion  is 
applicable 


biased”  beat  frequency  will  result  even  a  zero  rotation  rate  of 
the  laser  cavity. 

Let  6C w  and  0 Ccw  be  the  roundtrip  phase  shifts  of  the  CW 
and  CCW  waves,  L  be  the  roundtrip  optical  path  length,  and 
Xcw  and  Xccw  be  the  optical  wavelengths  in  the  laser  me¬ 
dium.  We  can  find  that  the  frequency  difference  of  the  CW 
and  CCW  waves  when  they  are  not  in  physical  rotation,  i.e., 
the  dc  bias  beat  frequency  due  to  nonreciprocal-diffraction 
filtering,  will  be  given  by 

( 2vL  n  \  /  2 vL  „  \ 

(r-  +<>cw  ‘  7-  +0CCw)-2H7r  (1) 

\Acw  /  \Accw  / 

or 

Aizs^cw  -  ('ccw)no  rotation 

=  £.  L  -  (gcw  ~  0ccw)1  ^ 

Lj  [_  2rr  J 

where  c  is  the  speed  of  light  inside  the  resonator,  and  n  is  the 
integer  difference  between  axial  mode  numbers  of  the  two 
waves.  The  lowest  beating  frequency  between  the  two  closest 
waves  is  obtained  by  choosing  the  proper  n  which  minimizes 
|Am|.  ' 

Wc  may  call  this  scheme  “diffraction  bias”  as  compared  to 
“Faraday  cell  bias”  or  “mechanical  bias.”  This  constant  dif¬ 
fraction  bias  may  be  able  to  get  rid  of  the  lock-in  phenomenon, 
in  exactly  the  same  way  as  a  Faraday  cell  docs  in  some  laser 
gyroscopes. 

This  diffraction  bias  has  the  following  properties. 

1)  An  intracavily  Faraday  cell  can  cause  problems  like  bias 
shifting,  etc.,  if  there  arc  magnetic  field  or  temperature  varia¬ 
tions.  Diffraction  bias  may  be  insensitive  to  changes  in  exter¬ 
nal  physical  parameters,  which  means  better  and  more  stable 
gyroscope  performance  can  l>c  achieved. 

2)  Compared  to  both  mechanical  dither  attd  Fatadav  cell 
bias,  dilTiaction  bias  may  be  simplei  in  conslioclion  and  per¬ 
haps  easier  in  maintenance. 

Astigmatism  isaneflccl  of  some  concern  in  ring  laser  design. 
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ccw 

OUTPUT 

Fig.  4.  With  designs  such  as  the  above,  keeping  ail  angles  of  incidence 
6  <  15°.  astigmatism  generally  will  be  small. 


(RADIUS  Rccw) 

Fig.  5.  Genera)  nomenclature  for  a  simple  negative  branch  unstable 
resonator. 


If  astigmatism  is  serious,  we  may  use  either  aspherical  mirrors 
or  a  small-angle  design  like  Fig.  4  to  minimize  its  effect.  Even 
with  astigmatism  our  former  discussion  is  still  true,  i.e.,  there 
are  still  nonoverlapping  focal  points  and  hence  differential- 
spatial  filtering  is  still  applicable  and  bias  can  be  obtained 
through  proper  implementation  of  filtering  pins.2 


II.  Geometrical  Analysis  of  Unstable  Ring 
Resonators 


We  will  first  develop  a  geometrical  analysis  to  predict  the 
locations  of  the  focal  points  in  an  unstable  ring  resonator.  Let 
Li,  L 2  be  the  separation  lengths  of  two  mirrors  with  radius  of 
curvature  Ft  and  F2,  respectively,  as  shown  in  Fig.  5.  For  the 
CW  wave,  we  choose  a  reference  plane  G  right  before  the  out¬ 
put  coupler  Rt.  The  roundtrip/lflCZ)  matrix  elements  for  the 
propagation  sequence  G-*F1->Lt->F2-+Z.2-*Gare  then  as 
follows: 


(3) 


At  plane  G  the  diverging,  perturbation-stable  geometric  eigen- 
modc  (2)  is  a  spherical  wave  with  radius  of  curvature  R</w 
given  by 


I 

/?(  w 


(A-D) 

/(*'*> Y  ,1 

2  \ 

'l  2  J  J 

(4) 


Ify?cw  >  o,  the  wave  is  divergent  in  the  direction  of  travel. 

The  CW  wave  is  focused  down  to  the  focal  point  Few  after 
reflection  off  the  output  mirror  /?,.  The  distance  Zcw  from 
R  i  to  Few  along  path  Lx  is  given  by 
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- + - = —  (5) 

ZcW  fll 


or 


Zcw 


*|/?cw 
2/?cw  "  Ri 


(6) 


assuming  that  0<ZCw  <Li,  so  that  Fcw  is  within  the  Lx 
segment.  If  few  falls  outside  this  range,  the  additional  focus¬ 
ing  effect  of  mirror  F2  must  be  taken  into  account.  (It  is  also 
possible  to  have  3,  5,  •  •  •  foci  in  a  more  complicated  multiele¬ 
ment  negative-branch  unstable  resonator,  but  here  we  are  only 
interested  in  the  single-focus  case.) 

For  the  CCW  wave  the  reference  plane  is  set  at  G'  and  all  the 
former  calculations  are  valid  with  the  exchange  of  L\  and  Z,2. 
The  CCW  focus  Fccw  is  now  located  at  /?ccw/2  behind  the 
output  mirror. 

A  typical  example,  with  confocal  design  in  the  CW  direction 
but  not  in  the  CCW  direction,  is  the  following: 


R\  =  18  cm 
F2  =  24  cm 
/,  i  =  21  cm 
1.2  ~  50.5  cm 

Z(-  w  =  4  cm  (on  /. ,  from  F, ) 


1  TTk-  jstipmatism  mentioned  here  i\  mainly  due  to  tilled  minors, 
where  the  cylindrical  symmetry  a  lonp  the  optical  path  is  spoiled  Asa 
result  we  have  2  two  told  symmetry  with  one  symmetry  avis  in  the 
plane  of  the  rinp  and  the  other  perpendicular  to  it.  Ibis  indicates  at 
least  two  pins  can  he  put  alone  either  axis  to  induce  nonreciprocal 
eipcnmodcs 


Zyi'w  =  b.l>  em  ("ii  /.  i  from  R , ) 

A  =  1.15  jr  (or  3.3*>#r) 

2d  =  2  mm  (oi  3.43  mm  for  3  3‘)  /a) 
(output  coupler  diameter) 
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7Veq  =  1 .5. 


normalized  radius 


Fig.  7.  (a)  Intensity  of  the  CW  I  -  0  eigenmode  versus  radius  (in  nor¬ 
malized  coordinates)  at  the  reference  plane  just  inside  the  output 
coupler,  (b)  Intensity  of  the  CCW  1=1  eigenmode  along  <>=0 
radius  at  the  other  reference  plane  just  inside  the  output  coupler. 


The  locations  of  Fcw ,  Fcc w  are  then  as  shown  in  Fig.  5.  The 
two  foci  are  separated  by  a  distance  of  2.1  cm. 

To  obtain  larger  focal  separation,  a  symmetric  confocal  de¬ 
sign  as  shown  in  Fig.  6  may  be  used.  At  these  focal  points, 
there  is  larger  mode  area  difference  between  the  CW  and  CCW 
waves  (see  following  analysis),  which  allows  easier  and  more 
efficient  nonreciprocal  filtering  to  exist. 

III.  Diffraction  Eigenmode  Analysis 


y&,(r)  =  exp  ^-/  (kL  -  j  •  exp  [-ftitDr*]  •  (tj)1 2rr 

r»lsf ififx 

•  I  fi/(fo)cxp  [-jtnArl] 

J 0 

•  Ji(2nrr0)  r„  dr0  (8) 

where  L  =  roundtrip  optical  length,  a  =  output  coupler  radius 


We  now  develop  a  full  diffraction  analysis  of  the  CW  and 
CCW  eigenmodcs  for  a  typical  example.  The  full  eigenmode 
integral  equation  for  the  unstable  resonator  with  roundtrip 
A,B,C,D  matrix  element  is3 


JJ  E(X°’K) 

'eXp[~rBx\{AxX°  '  **'*•.**■*'> 
exp  i  j-iAyY'o1  -  2Y'Y'0  +  DyY'1)] 


f 

r  o 

V\B\\ 


,  t  =  sgn  ( B ) 


and  /  =  0, 1 , 2,  3,  •  •  •  is  the  azimuthal  symmetry  number.  The 
extra  phase  shift  (tj)1  is  particularly  important  in  determining 
the  diffraction-bias  beat  frequency.  A  symmetric  confocal 
design  (Fig.  6)  was  tested  using  a  standard  resonator  program 
with  the  following  parameters: 


(9) 


■dx'ody'o.  (7) 

If  we  assume  cylindrical  symmetry  l’{X' ,  )  ')=  £/(/)  ■ 
and  change  to  llunkcl  transforms,  we  have 


I.\  =  /.j  =  41  cm 


3 The  propagation  phasor  is  assumed  In  be  r 
ill  signs  associated  with/  slmiild  he  changed 


If  ' 


is  given. 


R  |  =  35  cm 
Ri  =  47  cm 
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CCW  FOCAL  PLANE 
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2  3  4 

NORMALIZED  RAOIUS 


Fig.  8.  (a)  Intensity  of  the  CW  1=0  eigenmode  versus  radius  at  the  CW 
focal  plane  Few-  (b)  Intensity  of  the  CCW  /=  1  eigenmode  along 
C>  =  0  radius  at  the  same  Fcw  plane. 


Fig.  9.  (a)  Intensity  of  the  CW  1  =  0  eigenmode  versus  radius  at  the 
CCW  focal  plane  FqCV/-  (b)  Intensity  of  the  CCW  1  =  1  eigenmode 
along  0  =  0  radius  at  the  same  FqCVJ  plane. 


a  =  1  mm  (or  1 .72  mm  for  3.39  p) 

A=  1.15  n  (or  3.39  ft) 

2cw  =£ccw  =  1^-5  cm 

yv.q  =  -2.5. 

It  was  assumed  that  as  a  result  of  spatial  filtering  the  colli¬ 
mated  CW  wave  was  operating  in  the  /  =  0  mode,  and  the 
collimated  CCW  wave  was  filtered  and  operating  in  the  /  =  1 
mode.  An  iterative  computer  program  written  with  the  fast 
Hankel  transform  (FHT)  algorithm  (3]  was  used  to  find  the 
CW  and  CCW  eigenmodes,  which  are  shown  in  Fig.  7.  (Each 
mode  required  about  70  iterations  to  reach  the  convergence 
level  of  |  A7I  <  10~3.)  The  mode  patterns  at  the  CW  and  CCW 
foci  are  also  illustrated  in  Figs.  8  and  9.  These  clearly  show 
that  the  sizes  of  these  two  waves  are  significantly  different,  as 
expected.  Spatial  filtering  by  inserting  two  pins  to  constrain 
the  CCW  wave  but  not  the  other  should  definitely  be  possible. 
The  associated  complex  eigenvalues  in  the  two  directions  are 
found  to  be 

l7cwl(/-o)  *  0.888 
^'W(j.o)  =  1  5660 
l7ccwl(i-t)  =  0.846 

®l'cw</-i)  =  3.0013  +  — . 

The  two  oppositely  directed  modes  do  not  have  exactly  the 
same  losses  hut  the  losses  are  close  enough  that  simultaneous 
oscillation  in  the  two  directions  at  least  seems  conceivable, 
perhaps  with  further  adjustment  of  the  resonator  parameter 


The  extra  tr/2  of  0Ccw  «  due  to  the  phase  factor  (t/)'  in 
(8),  which  was  not  included  in  the  resonator  program.  The 
diffraction-biased  beat  frequency  for  this  case  is  Ap  1 76  MHz. 
Different  cavity  parameters  (e.g.,  /?, ,  R j,  L\,L j,  /)  were  also 
tried,  and  beat  frequencies  ranging  from  a  few  megahertz  to 
several  gigahertz  were  obtained.  More  systematic  information 
on  the  variations  of  transverse  beat  frequencies  in  unstable 
resonators  would  be  useful  in  this  connection. 

IV.  Conclusion 

The  idea  of  using  differential-spatial  filtering  or  “diffrac¬ 
tion  bias”  in  negative-branch  unstable  ring  resonators  to  obtain 
nonreciprocal  oscillation  frequencies  appears  to  be  a  new  idea. 
The  potential  application  of  this  technique  in  laser  gyroscopes 
is  particularly  interesting.  More  work  is  needed  to  see  what 
further  practical  advantages,  or  disadvantages,  are  associated 
with  this  technique. 
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We  outline  here  a  new  algorithm  for  evaluating  Hankel  < Four ler- Bessel)  transforms  numerically  with  enhanced 
speed,  accuracy,  and  efficiency.  A  nonlinear  change  of  variables  is  used  to  convert  the  one-sided  Hankel  transform 
integral  into  a  two-sided  cross-correlation  integral.  This  correlation  integral  is  then  evaluated  on  a  discrete  sam¬ 
pled  basis  using  fast  Fourier  transforms.  The  new  algorithm  offers  advantages  in  speed  and  substantial  advan¬ 
tages  in  storage  requirements  over  conventional  methods  for  evaluating  Hankel  transforms  with  large  numbers  of 
points. 


In  many  optical  calculations  one  would  like  to  evalu¬ 
ate  Hankel  transforms  numerically  with  computational 
efficiency  analogous  to  the  fast  Fourier  transform.  We 
describe  here  a  new  “fast  Hankel  transform”  (FHT) 
algorithm  that  can  evaluate  Hankel  transforms  while 
providing  many  of  the  desired  benefits. 

The  standard  Hankel  transform  of  order  /  is  given 
by 

gip)  =  2rr  rf(r)J,(2vpr)  dr  (1) 

with  a  symmetric  reverse  transform  integral  from  gip) 
to  fir).  Using  the  change  of  variables  r  =  r»e"x,  p  = 
Poe,,y,  where  r0,  p0,  and  o  arc  initially  arbitrary,  converts 
the  one-sided  Hankel  transform  integral  [Eq.  (1)]  into 
the  two-sided  cross-correlation  integral 

g(y)  =  f(x)j{x  +  y)dx,  (2) 

in  which  the  input  and  transform  functions  are  fix)  = 
r/(r),£(y)  =  pgip),  and  the  modified  Bessel  kernel  is  fix 
+  y)  =  Iirarp  Ji(2irrp).  This  change  of  variables  might 
be  called  the  “Gardner  transform.”  1  It  has  been  ap¬ 
plied  elsewhere  to  exponential  or  Laplace  transforms'2  1 
and  can  potentially  be  applied  to  any  integral  transform 
whose  kernel  has  a  product  argument/’ 

A  correlat  ion  integral  such  as  Eq.  2  can  be  evaluated 
efficiently  in  a  discrete  sampled  approximation  by  using 
Fourier-transform  methods.  Sampling  the  functions 
fix)  and  ffiy)  at  discrete  values  given  by  .t  =  n,  y  =  m 
withu.m  =0,1,...  ,7V  —  1  is  equivalent  to  sampling  the 
original  functions  fir)  and  gip)  at  discrete  exponentially 
increasing  values  given  by 

r„  =  roc"”,  p„,  =  pvf,,'m.  (3) 

These  sampled  values  are  confined  to  truncated  ranges 
r„  <  r  <  It  and  p,>  <  p  <  li  with  upper  limits  b  =  r„ r""v 
and  d  =  p„(",N.  The  product  yib  is  a  space  bandwidth 
product  for  the  transform.  In  terms  of  the  sampled 
values  f,(  -  fix,, )  =  r„fir„ ),  g,„  =  /?(ym )  =  p„,gip,„  I,  and 

. . =jixl;  +>’,„)  =  2k  <rr„p„,’J  ii'litr„pnl ),  the  Hankel 

transform  integral  |E<|.  (Dor  E(|.  (2)|  is  approximated 
to  some  level  of  accuracy  by  the  discrete  sum 


N—  1 

8m  ~  21  fnjn+m ■  '  (4) 

n  =  0 

If  the  input  sequence  /„  is  padded  with  zeros  and  the 
upper  limit  of  the  sum  extended  to  2 N  —  1,  this  sum 
becomes  a  27V-term  discrete  circular  correlation.  This 
discrete  correlation  can  be  evaluated  exactly  using  the 
discrete  2N-term  fast  Fourier  transforms  of  the  input 
sequences  f„  and  j„  by  carrying  out  the  steps 

gm  =  FFT  [FFT(/„, )  X  FFT*0* )|.  (5) 

The  FFT  and  FFT*  notations  indicate  the  forward  and 
inverse  (±j)  Fourier  transforms  (or  vice  versa).  The 
two  FFT's  inside  the  brackets  are  multiplied  together 
term  by  term  before  performing  the  outer  FFT. 

To  the  extent  that  the  values gm  generated  by  Eq.  (4) 
or  Eq.  (ft)  represent  a  good  approximation  to  the  true 
sampled  transform  values  ffiym ),  this  procedure  com¬ 
putes  an  N-point  discrete  approximation  to  a  Hankel 
transform  by  computing  two  2/V-term  fast  Fourier 
transforms  plus  2 N  additional  multiplications.  The 
transforms  are  performed  in  place,  requiring  only  the 
additional  storage  of  2 TV  previously  computed  values 
of  the  Fourier  transform  of The  extra  values  of 
gm  from  m  =  N  to  m  =  2 N  —  1  represent  aliased  results 
from  the  circular  correlation  process  and  are  discard¬ 
ed. 

Criteria  for  choosing  the  parameters  r)t,  po.  and  «v  can 
be  developed  as  follows.  If  d  is  the  largest  value  of  p  in 
the  transform  domain,  the  highest  frequency  compo¬ 
nent  in  the  r  domain  is .7,(27 rdr)  ~  cos  (2;rdr).  We  then 
require  that  the  lower-end  truncation  point  in  the  r 
domain  not  he  larger  than  MK\  cycles  of  the  highest 
spatial  frequency  d.  or  r()  =  (K i/7)“'.  Similarly  we  re¬ 
quire  that  the  sample  point  spacing  Ar„  '  r„+ ,  -  r„  at 
the  upper  truncation  point  r  =  h  not  be  greater  than 
\/K<  cycles  of  d.  or  Sr\  -  o/>  =  il\  >ii)'  ’.  Combining 
these  arguments  with  exactly  analogous  arguments  in 
the  p  domain,  namely  p,t  ~  il\ih)~]  and  Sp\  =  o d  = 
(/\ _./>)'  '.  leads  to  the  set  of  relations 

N  =  K-jiih  ln(/\,d/>), 
ar"N  =  l\,/H>, 

r  ii/<o  =  il\/K  I'-’lir.  ((>) 
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(a)  (b) 

Fig.  1  fa)  Exact  analytical  functions  fix),  j(x),  and  g{x) 
corresponding  to  a  uniform  circular  “top-hat”  input  function 
fir)  and  its  Airy -disk  Hanket  transform,  including  the  top-hat 
input  function  /  (top  curve);  the  /  =  0  Bessel-function  kernel 
j  (middle  curve);  and  the  Airy-disk  Hankel  transform  g 
(bottom  curve).  Parameter  values  for  the  change  of  variables 
are  a  =  0.01612,  r,,  =  />()  =  0.06349,  and  ft  =  b  =  3.938,  corre¬ 
sponding  to  N  =  256,  K  i  =  K  i  =  4.  The  "dimples"  in  the  g 
function  at  larger  arguments  are  spurious  results  of  a  plotting 
routine  with  inadequate  point  spacing,  (b)  Discrete  Hankel 
transformation  and  then  backtransformation  of  a  sampled 
top-hat  input  sequence  /„  (top  curve)  into  the  transformed 
sequence  g„  (middle  curve)  and  then  back  into  a  sequence  fn 
(bottom  curve),  using  the  FHT  algorithm  with  N  =  256,  K\ 
-  Ki-  4.  The  portions  of  the  transformed  and  backtrans- 
formed  sequences g„  and  f„  above  n  =  256  arise  from  aliasing 
effects  in  the  transform  algorithm  and  are  discarded.  The 
backtransformed  sequence  f„  exhibits  Gibbs  phenomena 
because  of  the  finite  truncation  of  theg„  sequence  before  it 
is  backtransformed. 


We  hypothesize  that  the  "points-per-cycle”  parameters 
K |  and  K ■>  should  both  be  2:2  points  per  cycle.  Equa¬ 
tions  (6)  then  determine  a,  r(1/>n,  and  fib  appropriate  to 
a  given  number  of  points  N;  or  alternatively  the  number 
of  points  N  and  the  a  and  r(>p<>  parameters  required  to 
transform  with  adequate  accuracy  a  function-transform 
pair  having  a  given  space-bandwidth  product  (ib. 

This  algorithm  was  programmed  in  KOliTRAN  using 
single-precision  arithmetic  on  an  IBM  370/168  com¬ 
puter.  Figure  1(a)  shows  for  purposes  of  comparison 
analytical  plots  of  a  "top-hat"  input  function /(r)  =  1, 
0  <  r  <  b,  transformed  into  /(*):  the  modified  Bessel 
kernel j(x);  and  (he  Airy-disk  function/?!*)  that  is  the 
exact  Hankel  transform  of  a  uniform  input;  all  plotted 
as  continuous  functions  for  parameter  values  corre¬ 
sponding  to  I  =  0.  .V  =  256,  and  K i  =  K  >  =  4  points/ 
cycle.  The  kernel  function  jix)  is  seen  to  be  slowly  di¬ 
vergent  at  large  ,r. 

Figure  1(b)  shows  the  discrete  numerical  results 
generated  by  the  FHT  algorithm  for  this  case.  The 
FI  IT  algorithm  samples  the  input  fund  ion  fix )  |or  at 
lernalivelv  .!?( x  (|  at  integer  values  from  I)  to  255.  pads 
the  input  sequence  with  zeros  from  256  to  5|  I.  performs 
a  discrete  circular  convolution  against  512  sampled 
points  from  the  kernel  /I  \  ).  and  yields  an  output  se¬ 
quence  that  approximates  the  exact  Hankel  transform 


over  the  range  0-255  and  contains  aliased  results  over 
the  range  256  51 1. 

In  Fig.  1(b)  the  middle  curve  shows  the  output  se¬ 
quence  f>m  generated  by  the  FHT  algorithm  for  the 
“top-hat”  input  sequence  of  the  top  curve,  using  N  = 
256  and  K\  =  K->  —  4.  The  Airy  disk  transform  gift)  = 
2Ji(2irb(i)/i>  is  accurately  reproduced  by  the  FHT  al¬ 
gorithm  for  0  <  m  <  N,  with  irrelevant  aliased  values 
for  m  >  /V.  The  sequence  g,„,  when  truncated  and  fed 
back  through  the  FHT  algorithm,  provides  a  lair  re¬ 
production  of  the  original  top  hat,  as  shown  in  t  he  bot¬ 
tom  curve;  but  with  obvious  Gibbs-phenomena  ripples 
resulting  from  the  finite  truncation  of  the  Airv-disk 
transform. 

The  Laguerre-Gaussian  functions  upiir)  =  (2irr)- 
Lpl(2irr-)  exp(— irr-)  are  self-transforming  under 
Hankel  transformation.  Figure  2  shows  in  r.p  coordi¬ 
nates  the  FHT  algorithm  applied  twice  in  succession  to 
a  Laguerre-Gaussian  function  with  /  =  0  and  p  =  8, 
using  N  =  128  with  K\  =  K>  =  2.  The  mean-square 
error  on  successive  transforms  in  this  case  is  5  0.4%. 
Laguerre-Gaussians  have  been  accurately  transformed 
out  top*  100  for  N  =  1024,  K  ]  =  8,  and  K±  =  2. 

Figure  3  shows  a  top  hat  and  also  an  annular  function 
transformed  and  then  backtransformed  using  the  FHT 
algorithm.  The  observed  Gibbs-phenomena  ripples  are 
expected  for  finite  truncation  in  the  transform  domain 
whatever  Hankel  transform  algorithm  may  be  em¬ 
ployed.  They  can  presumably  be  reduced  by  using 
proper  windowing  (apodizing).  The  "droop"  in  the 
function  values  at  low’  r  is  understood  and  can  be  elim¬ 
inated  also. 

A  small  number  of  experiments  indicate  that  the 
FHT  algorithm  works  equally  well  for  higher-azi¬ 
muthal-order  Hankel  transforms.  Beyond  the  pre¬ 
liminary  tests  summarized  here,  the  capabilities  of  this 


Fig.  2  The  input  function  and  two  successive  KMT's  of  a 
Daguerre  (laussian  input  function  with  p  =  M,  I  =  n  using  the 
Fit T algorithm  with  ,V  =  IJN.  A,  =  A  .  =  2.  The  lirst  FIIT 
lies  esse n I  iallv  on  lop  ot  t  he  input  fund  ion  at  low  arguments 
and  I  hen  drops  slight  l\  lower.  The  second  FHT  i-  uni  iallv 
slightly  lower,  then  follows  closely  the  first  Fill'  the  curves 
are  plotted  as  straight  lines  between  the  discrete  sample 
points;  the  discrete  |mml  spacing  is  evident  at  larger  argu 
incuts,  s 


,1 


-  45 


July  1977  /  Vol.  1.  No.  1  /  OPTICS  LETTERS  .  15 


Fig.  3  Results  of  Hankel-transforming  and  then  back- 
transforming  a  top-hat  input  function /(r)  =  1,0  <  r  £  6,  and 
a  uniform  annular  function  fir)  =  1, 6/4  <  r  <  36/4,  using  the 
FHT  algorithm  with  N  =  512,  K\  =  K_>  =  4  («  =  0.0091648,  r <> 
=  0.0478665,  6  =  5.223).  The  slump  in  the  top-hat  response 
and,  to  a  lesser  extent,  in  the  annular  response  at  low  r  is  re¬ 
lated  to  the  truncation  at  p0  of  the  narrow  Airy-disk  Hankel 
transform  in  the  p  domain;  rescaling  of  the  r  and  p  coordinates 
to  give  comparable  widths  in  the  two  domains  gives  much 
better  results. 


algorithm  remain  to  he  fully  determined  (and  opti¬ 
mized).  One  can  compare  the  computational  efficiency 
of  the  FHT  algorithm  with  straightforward  numerical 
evaluation  of  Hankel  transforms  using  trapezoidal  in¬ 
tegration  at  uniformly  spaced  sample  points  with  K 
points  per  cycle  in  both  r  and  p  domains.  Assuming 
that  the  Bessel  function  values  will  he  stored  in  a  look¬ 
up  table,  this  calculation  will  require  N  =  Kdb  sample 
points  in  each  domain  and  VjA'’-  =  '/>(Kd6)2  complex 
multiplications  per  transform.  The  FHT  algorithm, 
including  the  aliased  values,  requires  a  larger  number 
2  N  =  2  K>fib  ln(/f  i//6)  of  sample  points  because  of  the 
unavoidable  oversampling  with  exponentially  spaced 
sample  points.  However,  only  ~2N  log_>  4/V  complex 
multiplications  are  required  to  compute  the  FFT  s  plus 
additional  steps  in  the  FHT  algorithm.  The  FHT  al¬ 
gorithm  requires  significantly  fewer  operations  for 
space- bandwidth  products  above  fib  '  140  for  K  =  K\ 
=  Kj  =  2  or  above  fib  *  70  for  h  =  l\  \  =  K  >  =  4  points 


per  cycle.  One  might  also  hope  to  obtain  improved 
accuracy  against,  round-off  errors  because  of  the  supe¬ 
rior  properties  of  the  FFT  algorithm,  alt  hough  t  his  point 
has  yet  to  be  tested. 

Of  greatest  practical  importance,  however,  is  the  re¬ 
duction  in  computer  storage  requirements  using  the 
FHT  algorithm.  Straightforward  Hankel-transform 
evaluation  using  stored  Bessel-lunction  values  requires 
storing  lf2(Ktib)-  real  values  of  the  Bessel  kernel 
Ji(2irrp)  for  equally  spaced  values  of  both  r  and  /<  (or 
else  the  use  of  interpolation  techniques  that  require 
additional  multiplications).  The  FHT  algorithm,  in 
contrast,  is  carried  out  in  place,  requiring  storing  only 
8K>f1b  ln(A',//6)  total  real  numbers  in  the  function  and 
kernel  arrays.  For  fib  =  1000  and  K  ]  =  K  >  =  2,  this 
would  require  ~500  Kbytes  of  storage,  assuming  4  bytes 
per  real  number,  as  against  ~8000  Kbytes  in  the  stan¬ 
dard  case. 

These  advantages  in  speed,  storage  requirements,  and 
potential  accuracy  should  make  the  FHT  algorithm 
valuable  in  optical-beam  and  resonator  calculations  as 
well  as  in  other  optical  and  nonoptical  applications. 
Further  tests  and  improvements  on  the  basic  algorithm 
will  be  reported  in  future  publications. 

This  work  was  carried  out  with  support  from  KR1)A 
through  the  Los  Alamos  Scientific  Laboratory;  from 
United  Technology,  Pratt  &  Whitney  Aircraft;  and  from 
the  Air  Force  Office  of  Scientific  Research. 
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ORTHOGONALITY  PROPERTIES  OF  OPTICAL  RESONATOR  EIGENMDDES 
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ABSTRACT 


The  transverse  eigenmodes  of  optical  resonators  are  not  in  general 
orthogonal  to  each  other,  but  rather  are  biorthogonal  to  an  adjoint  set 
of  functions  which  turn  out  to  be  the  transverse  eigenmodes  for  propa¬ 
gation  in  the  reverse  direction  around  the  same  resonator.  We  present 
a  fairly  general  derivation  of  this  property  which  is  valid  for  ring 
or  standing-wave  and  for  stable  or  unstable  resonators  with  arbitrary 
output  mirror  profiles.  We  also  point  out  certain  unusual  consequences 
of  these  properties  for  high -magnification  unstable  resonators. 
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ORTHOGONALITY  PROPERTIES  OF  OPTICAL  RESONATOR  EIGENMODES 

by 

A.  E.  Slegman 

Edward  L.  Ginzton  Laboratory 
Stanford  University 
Stanford,  California  USA 

I .  INTRODUCTION 

The  orthogonality  properties  of  optical  resonator  transverse  modes  ap¬ 
pear  to  be,  if  not  inadequately  understood,  at  least  not  very  widely  under¬ 
stood  and  very  little  discussed  in  the  literature  [1-3] .  The  basic  operators 
that  determine  optical  resonator  transverse  eigenmodes  are  in  general  not 
hermitian,  and  therefore  these  eigenmodes  are  not  in  general  orthogonal  in 
the  (u*, u)  sense,  where  u  means  complex  conjugate.  Consider  first  a 
simple  two-mirror  standing-wave  resonator  with  arbitrarily  shaped  but  uni¬ 
formly  reflecting  mirror  surfaces.  The  Huygens  kernel  for  the  standard 
"Fox-and-Li"  approach  to  finding  the  resonator  eigenmodes  is  then  symmetric 
between  points  on  the  mirror  surfaces.  The  resonator  eigenmodes  in  this  case 
are  therefore  orthogonal  in  the  (u, u)  sense  over  the  mirror  surfaces  [J] . 
The  fundamental  adjoint  properties  of  the  wave  equation  and  its  associated 
Green's  function  [4]  strongly  indicate  that  this  orthogonality  property 
probably  applies  much  more  generally  to  a  broader  class  of  resonators. 
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Arnaud  [5]  discusses  the  biorthogonal  eigenmodes  of  general  complex- 
paraxial  resonators^  stable  or  unstable;  having  general  hermite-gaussian 
modes.  He  points  out  that  there  is  a  set  of  adjoint  eigenmodes  u’  which 
are  biorthogonal  to  the  usual  set  of  resonator  eigenmodes  u  in  this  case; 
and  that  the  adjoint  modes  are  physically  the  transverse  eigenmodes  for 
propagation  in  the  reverse  direction  around  the  same  resonator.  Again  it 
seems  likely  that  this  result  applies  more  generally  to  nearly  all  optical 
resonators. 

The  purpose  of  this  note  is  to  present  a  simple  derivation  which  demon¬ 
strates  the  biorthogonality  properties  of  the  eigenmodes  of  an  optical  reso¬ 
nator  under  still  more  general  conditions;  including  an  arbitrary  output 
mirror  profile  in  both  phase  and  amplitude.  The  conclusions  obtained  then 
apply  to  stable  and  unstable  resonators,  to  ring  and  standing-wave  resonators, 
and  especially  to  resonators  with  non-paraxial  and  non-uniform  elements 
located  at  least  one  plane  within  the  cavity.  We  also  note  that  these 
biorthogonality  properties  imply  certain  unusual  consequences  for  high- 
magnification  unstable  resonators. 
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RESONANT  MODES  OF  OPTICAL  CAVITIES 
WITH  PHASE-CONJUGATE  MIRRORS 

Pierre  A.  Be'langer,  Amos  Hardy,  and  A.  E.  Siegman 

Edward  L.  Ginzton  Laboratory 
Department  of  Electrical  Engineering 
Stanford  University  . 

Stanford,  California 

(Revised  25  October  1979) 

ABSTRACT 

We  derive  the  lowest-order  self-consistent  gaussian  transverse  modes  and 
also  the  resonant  frequencies  of  an  optical  resonator  formed  by  conventional 
paraxial  optical  components  plus  a  phase-conjugate  mirror  (PCM)  on  one  end. 
The  conventional  optical  elements  are  described  by  an  overall  ABCD  matrix. 
Cavities  with  purely  real  elements  (no  aperturing)  have  a  continuous  set  of 
self-reproducing  gaussian  modes,  described  by  a  semicircular  locus  in  the 
l/q  plane  for  one  round  trip;  all  gaussian  beams  are  self-reproducing 
after  two  round  trips.  Complex  ABCD  matrices,  such  as  are  produced 
by  gaussian  aperturing  in  the  cavity,  lead  to  unique  self-consistent 
perturbation-stable  gaussian  modes.  The  resonant  frequency  spectrum 
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of  a  PCM  cavity  consists  of  a  central  resonance  at  the  driving  frequency 

of  the  PCM  element,  independent  of  the  cavity  length  L  ,  plus  "half- 
axial"  sidebands  spaced  by  &oax  *  27t(c/^l)  ,  with  phase  and  amplitude 
constraints  on  each  pair  of  upper  and  lower  sidebands. 
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I.  INTRODUCTION 


The  emerging  subject  of  phase  conjugate  optics  is  receiving  much  attention 

12  3 

at  present.  *  * ^  The  basic  building  block  for  phase  conjugate  optics  is  the 
phase  conjugate  mirror  (PCM),  an  optical  element  that  has  the  unusual  property 
of  converting  an  incident  optical  wave  into  a  reflected  wave  whose  complex 
phasor  amplitude  is  the  complex  conjugate  of  the  incident  wave.  In  practice 
phase-conjugate  reflection  is  most  often  accomplished  by  means  of  four-wave 
mixing  in  a  nonlinear  optical  medium,  such  as  an  optical  Kerr  cell,  when 
this  medium  is  actively  pumped  by  two  coherently  related,  counter-propagating 
pump  beams.  Phase-conjugate  reflection  can  also  be  accomplished  using  three- 
wave  mixing  (degenerate  parametric  mixing)  but  with  severe  practical  limi¬ 
tations  set  by  phase-matching  requirements. 

Since  the  phase  conjugate  mirror  is  an  interesting  and  unusual  new  com¬ 
ponent,  our  objective  in  this  paper  is  to  explore  the  transverse  and  longi¬ 
tudinal  resonant  mode  properties  of  optical  cavities  having  phase  conjugate 

mirrors  as  one  of  their  end  reflectors.  To  our  knowledge  only  one  laser  has 

4 

yet  been  operated  using  such  a  pure  PCM  end  mirror.  Several  other  analyses 

4. 5. 6 

of  this  topic  have  appeared  while  this  paper  was  in  preparation.  ’  Our 
results  here  partially  overlap  and  partially  complement  and  extend  these 
contributions . 

In  this  paper  we  first  carry  out  in  Section  II  a  general  derivation  of 
the  lowest-ordcr  gaussian  transverse  mode  properties  of  a  PCM  laser  cavity. 


-  53  - 


r 

We  include  the  possibility  of  radially  varying  transmission  or  gain,  which 
leads  to  complex-valued  paraxial  elements.  Section  III  then  explores  a 
few  specific  examples  of  the  transverse  modes  in  simple  PCM  cavities  con¬ 
taining  either  gaussian  apertures  (i.e.,  thin  elements  with  gaussian  radial 
apodization)  or  gaussian  ducts  (i.e.,  thick  elements  with  gaussianly  varying 
transmission).  Finally  in  Section  IV  we  discuss  the  somewhat  unusual  resonant- 
frequency  or  axial-mode  characteristics,  particularly  the  existence  of  "half- 
axial"  modes,  that  are  characteristic  of  PCM  laser  cavities. 
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II.  GAUSSIAN  TRANSVERSE  MODE  ANALYSIS 


A.  Round-Trip  Gaussian  Beam  Propagation 

We  consider  as  our  basic  model  an  ideal  PCM  reflector  placed  at  one  end 
of  a  general  multi-element  optical  resonator,  as  in  Fig.  1.  A  PCM  reflector 
has  the  basic  property  that  if  the  optical  wave  incident  on  the  reflector 
along  the  positive  z  axis  has  the  form  at  any  reference  plane  z 

~  J  (u>0t  -  kz  ) 

einc(x,y,  z,t)  =  Re  E(x,  y,  t)  e  (l) 

\ 

where  x,y  are  transverse  coordinates  across  a  reference  plane,  then  the 
reflected  wave  at  the  same  plane  has  the  form 

^  j  (coQt  +  kz) 

Cre{i(x,y,z,  t)  ■  Re  K  E  (x,y,t)  e  (2) 

where  k  is  a  PCM  reflection  coefficient  which  may  be  complex  (i.e.,  have 
a  phase  angle).  The  frequency  u^j  is  the  reference  optical  frequency  given 
by  the  frequency  of  the  pump  lasers  necessary  to  drive  the  nonlinear  medium 
used  in  the  PCM. 

As  a  consequence  of  the  phase  conjugation,  a  gaussian  optical  beam  bounc¬ 
ing  off  a  rcM  will  have  its  radius  of  curvature  R  reversed  in  sign,  but  its 
spot  size  w  unchanged  (except  for  a  potential  generalization  to  be  discussed 
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later).  If  we  define  the  complex  q  parameter  for  a  gaussian  beam  in  the 
usual  way,  i.e. 


then  the  reflective  property  of  the  PCM  in  Fig.  1  can  be  expressed  mathematically 
by 


The  other  optical  elements  in  the  resonator  outside  of  the  PCM  can  be  de¬ 
scribed  by  a  general  paraxial,  possibly  complex-valued  ABCD  matrix  which  sums 
up  a  complete  round-trip  from  the  reference  plane  out  through  the  resonator 
and  back  to  the  PCM.  Since  this  is  an  "out-and-back"  journey,  the  ABCD 
matrix  has  the  general  property  that  A  =  D  .  We  assume  negligible  dif¬ 
fraction  effects  at  any  apertures  or  stops  inside  the  resonator. 

A  gaussian  beam  with  a  general  complex  gaussian  beam  parameter  when 

it  leaves  the  PCM  will  then  return  with  a  gaussian  parameter  q^  given  by 

Aqg  +  B 

q'  „  - 

Cq2  +  D 

Therefore  a  complete  round  trip  around  the  entire  FCM  cavity  may  be  expressed 
by  the  net  transformation 


B  A2  -  1  -  AB/q* 

t3  — . ,  —  , .7. 

qj  A  -  B/q* 


(0) 
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where  we  have  used  the  fact  that  the  paraxial  ray  matrix  will  obey  the  re¬ 
lationship  AD  -  BC  *  1  • 

The  self-consistent  or  self-reproducing  solutions  to  Eq.  (6),  defined 
by  the  condition  that  **  ■  q  ,  are  then  given  by  the  solutions  to 


-t  +  4)  +  a2-  i  -  o 

qq  \q  q  / 


(7) 


This  defines  either  a  circular  locus  or  a  unique  singular  point  in  the  l/q 
plane,  depending  only  on  whether  the  ABCD  coefficients  are  purely  real 
or  complex. 


B.  Real  Paraxial  Systems 


For  purely  real  ABCD  elements  Eq.  (7)  is  a  purely  real  equation,  which 
may  be  recast  into  the  useful  form 


(8) 


This  defines  a  displaced  unit  circle  in  the  complex  B/q  plane,  that  is, 
in  a  plane  with  rectangular  coordinate  axes  x  s  b/r  and  y  £  Bx/nv2  , 
as  shown  in  Fig.  2(a).  This  displaced  unit  circle  has  radius  unity  and 
center  at  b/r  =  A  ,  Bx/nw2  “0  .  There  is  evidently  no  unique  self- 

consistent  transverse  mode  for  the  PCM  cavity  with  purely  real  elements. 
Any  gaussian  beam  corresponding  to  a  point  on  the  semicircular  locus 
represents  an  equally  valid  self-consistent  confined  mode. 
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Physical  insight  into  this  situation  can  be  obtained  by  examining  the 
variation  of  an  arbitrary  gaussian  beam  on  two  successive  round  trips,  start¬ 


ing  from  any  arbitrary  q  value  determined  by  an  initial  R  and  w  , 
located  for  example  inside  the  self-consistent  locus.  The  complex  vector 
(B/q  -  a)  from  the  center  of  the  circle  (a,  0)  to  the  initial  point  B/q 
in  the  complex  plane,  and  the  corresponding  vector  (ll/q'-A)  after  one 
round  trip  are  then  related  by 


(B/q  -  A) (B/q'  -  A)*  ■  1 


(9) 


This  describes  geometrically  an  inversion  with  respect  to  the  displaced 
unit  circle  in  the  complex  B/q  plane,  as  illustrated  in  Fig.  2(b).  This 
necessarily  implies  that  a  second  round  trip  will  reinvert  the  q'  parameter 
back  to  its  original  q  starting  place.  Hence,  with  real  ABCD  elements  all 
gaussian  beams  are  self-reproducing  after  two  round  trips.  The  parameter  q 
(or  B/q)  bounces  back  and  forth  between  conjugate  points  inside  and  outside 
the  self-consistent  locus,  as  illustrated  in  Fig.  2(b). 

This  behavior  occurs  because  the  PCM  acts  in  essence  as  an  automatically 
adjusting  lens  which  exactly  refocuses  a  wave  so  as  to  simply  change  the 
sign  of  R  .  Figure  J  shows  an  unfolded  model  for  a  PCM  cavity,  illus¬ 
trating  the  automatic  self-consistency  after  two  round  trips. 

We  may  also  conclude  that  the  confined  self-consistent  solution  at  any 
point  on  the  displaced  semicircle  is  marginally  stable  ~  which  means  it  is 

effectively  stable  —  against  perturbations  in  the  sense  introduced  by 
7 

Caspcrson.  This  occurs  because  any  input  beam  that  is  displaced  very 
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slightly  from  any  point  on  the  self-consistent  locus  will  bounce  back 
and  forth  across  the  locus,  without  any  growth  of  the  displacement. 

The  type  of  self-consistent  mode  analysis  we  have  carried  out  here 

using  a  reference  plane  in  front  of  the  PGM  can,  alternatively,  be  carried 

out  using  a  reference  plane  located  on  the  opposite  end  mirror  surface  of 

the  cavity,  e.g.  on  the  left-hand  mirror  surface  in  a  cavity  like  Fig.  1. 

The  result  of  such  an  analysis  shows  that,  for  real  ABCD  elements,  the 

phase  front  of  the  gaussian  beam  exactly  matches  the  mirror  surface  (the 

wave  radius  R  coincides  with  the  mirror  radius  of  curvature);  while  the 

spot  size  w  on  the  mirror  surface  is  arbitrary.  This  property  of  PCM 

4 

resonator  modes  was  pointed  out  earlier  by  AuYeung.  The  locus  of  q 
values,  referenced  to  the  end  mirror  surface,  is  thus  a  vertical  line 
in  the  l/q  plane  of  Fig.  2.  The  one-way  propagation  from  this  mirror 
through  the  optical  elements  to  the  PCM  end  of  the  cavity  is  a  biorthogonal 
transformation  in  the  l/q  plane  which  converts  this  vertical  line  into  the 
circular  locus  of  Fig.  2  as  the  spot  size  on  the  end  mirror  surface  varies. 

The  allowed  values  of  the  gaussian-beam  q  parameter  are  thus  dis¬ 
tributed  along  a  continuous  locus  for  real  ABCD  values,  rather  than 
being  limited  to  a  single  point.  However,  it  seems  very  likely  both  from 
conclusions  we  will  derive  in  the  following  section  and  from  the  analysis 
of  Bel'dyugin,  Galushkin,  and  Zanskov, ^  that  if  any  kind  of  weak  aperture 
is  added  inside  the  PCM  cavity,  diffraction  effects  will  remove  this  de¬ 
generacy.  The  cavity  will  then  converge  to  a  unique  quasi -gaussian  mode 
near  that  point  on  the  semicircular  locus  which  minimizes  the  diffraction 
losses  at  the  aperture. 
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Complex  Paraxial  Systems 


C. 

In  the  most  general  paraxial  case  an  optical  cavity  may  Include  among 
its  elements  a  "gaussian  aperture",  by  which  we  mean  a  thin  apodlzed  element 
with  a  radial  amplitude  transmission  varying  as  exp  (-  QSc^)  ,  and/or  a 
"gausslan  duct",  by  which  we  mean  a  thick  element  with  the  same  property. 

Such  elements  can  be  described  as  complex  paraxial  elements  —  that  is,  all 
of  the  rules  of  paraxial  optics  with  ABCD  matrices  remain  valid,  but  at 
least  some  and  potentially  all  of  the  ABCD  matrix  elements  become  complex¬ 
valued  quantities. 

It  becomes  convenient  in  this  case  to  write  the  ABCD  matrix  elements  in 
the  general  form 


A 

B 


B0e 


(10) 


where  Aq  and  Bq  are  the  absolute  magnitudes  of  the  complex  elements  A 

and  B  .  ■  The  angles  6  and  0_  in  general  have  finite  values  determined 

A  B 

by  the  gaussian  variable  amplitude  transmission  elements  inside  the  optical 
system.  These  angles  go  either  to  0  or  to  7T  as  the  strength  of  the 
gaussian  aperturing  goes  to  zero  (the  value  tt  meaning  that  the  matrix 
element  becomes  real  but  negative). 


We  note  here  that  a  gaussian  transverse  variation  of  the  amplitude  trans¬ 
mission  may  be  produced  by  the  PCM  mirror  itself.  That  is,  the  4-wave  mixing 
cell  used  to  produce  the  PCM  may  be  pumped  by  pump  beams  which  have  gaussian 
transverse  amplitude  variations;  and  this  in  turn  will  produce  a  gaussian 
transverse  variation  in  the  magnitude  of  the  PCM  reflection  coefficient 

K  .  The  net  effect  will  be  the  same  as  if  a  gaussian  aperture  were  added 

q 

to  the  optical  system  just  in  front  of  the  PCM.  ° 

The  self-consistency  equation  (7)  now  becomes  complex,  with  finite  real  and 
imaginary  parts.  These  yield,  after  a  little  algebra,  the  unique  self-consistent 
solutions 


B  sin  2 0 

— 2-  *=  A-  cos  (0  -  0  J - 

Rsc  ^0  Sin(0B  - 0A) 


(U) 


and 


m  - 

\7IW  / 

'  sc  ' 


2  r  .2  ,2 


[Ag  Sin*  (0B  -  0A)  -  Sin2  0B]  [Aq  sin2  (0fi  -  0A)  +  cos*  0^] 

4  Sl“2  1%  -  V 


(12) 


22  1  sin  20 

cos  20g  +  Aq  sin  (0fi  -  ©A)  1 - 


■( 


2A0  sin  (0g  -  0A) 
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These  results  say  that  if  the  condition  Is  satisfied  that 


sin2  (eg  -  eA)  2 


2 

sin  0 
_ I 

4 


(13) 


2 

then  a  complex-valued  paraxial  system  will  have  a  unique,  confined  (w  >  0), 
self-consistent,  lowest-order  gaussian  mode,  corresponding  to  a  point,  not  a 
curve,  in  the  B/q  plane.  As  0^  and  0^  approach  zero,  this  point  ap¬ 
proaches  in  the  limit  a  location  somewhere  on  the  semicircular  locus  for  real 
ABCD  systems.  However,  just  where  this  limiting  point  falls  on  the  displaced 
unit  circle  locus  depends  on  Aq  and  on  the  relative  values  of  the  angles 

0.  and  0„  ,  as  we  will  illustrate  in  the  following  section. 

A  B 

Suppose  in  particular  that  0  and  0  are  both  small,  approaching  zero 

A  "  B 

in  the  limit.  Then  the  self-consistent  solution  is 


B 

“  Aq 

V 

Rsc 

WV 

(1°L\ 

-  1  -  1 

f  ) 

-  v  / 

(I*) 


To  this  level  of  approximation  the  complex  self-consistent  point  lies  on  the 
displaced  unit  circle,  at  an  angle  cp  from  the  x  axis  given  by 


cos  <p 


V°b  -  V 


(15) 


as  illustrated  in  Fig.  H. 


We  must  also  check  whether  the  complex  self-consistent  solution  given  by 

*T 

( 11),  and  (12)  is  perturbation-stable  in  the  sense  described  by  Casperson.  1 

If  a  gaussian  beam  starts  around  the  PCM  cavity  with  an  initial  gaussian 

beam  parameter  given  by  q  +  Aq  ,  where  q  is  the  self-consistent  value 

given  by  (ll)  and  (12),  and  Aq  is  a  small  complex  perturbation,  then  after 

one  round  trip  the  beam  will  be  described  by  q  +  Aq'  where  Aq'  is  a 

s  c 

slightly  different  perturbation.  Perturbation  stability  requires  that 
|Aq'/Aq|  £  1  •  For  the  PCM  resonator  it  can  be  shown  that  after  one 

round  trip 

Aq'  2  [Re(A/B  -  l/q  )]2  +  [I©(a/b  -  l/q  )]2 
_ _  s>  c _ s  c 

Aq  [Re(A/B  -  l/qsc)]2  +  [I*U/B  +  l/qsc)]2 

[Aq  cos  (eB-0A)  ~  Bq/r sc]2  +  [A0  sin  (0B-0A)  -  Bqx/tiw2c]2 

(Aq  cos  (0b  -  eA)  -  Bq/Rsc]2  +  [AQ  sin  -  0^  +  BqX/tiw^]2  (l6) 

This  is  less  than  or  equal  to  unity  if 


Im  (a/b)  <;  0 

(17) 

sin  (0  -  e  )  £  0 

B  A 

(18) 

Figure  4  illustrates  how  a  stable  system,  starting  from  an  arbitrary  initial 
q  value,  will  converge  in  to  the  self-consistent  value  q  upon  repeated 
bounces . 
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III.  EXAMPLES  OF  COMPLEX  PCM  RESONATORS 


A.  Plane-Mirror  Cavity  With  Gaussian  Aperture 

As  a  particular  example  of  a  complex  PCM  resonator  let  us  consider  a 

cavity  consisting  of  a  plane  mirror  and  a  PCM  separated  by  distance  L  ,  plus 

a  gaussian  aperture  with  amplitude  transmission  T(r)  =  exp  (-  r  /a  )  located 

at  a  distance  L  from  the  plane  mirror,  as  in  Fig.  5(a) •  We  may  define  an 
a 

aperture  strength  parameter  F  and  a  normalized  aperture  location  X  by 

■na2  “  L 

F  s  -  ,  X  =  -S-  (19) 

LX  L 

The  complex  paraxial  matrix  elements  are  then  given  by 


A  »  1  -  2X(l  -  X)/F2  -  j  2/f 
B  =  2L  [1  -  X(l  -  X)2/F2  -  j  (1~X2)/f] 


(20) 


In  the  limit  of  weak  aperturing  (F  ->«)  this  gives 


Awl  ,  B  «  2L 

o  ’  o 

0a  w  -  2/F  ,  «  -  (1-X2)/F 


(21) 
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The  self-consistent  solutions  at  a  reference  plane  just  Inside  the  PCM  are 


sc 

LX 


1  +  X 
X 


(22) 


p  2 

7TW  1  +  X 

sc 

These  solutions  move  along  a  contour  in  the  l/q  plane  as  in  Fig-  5(h). 

As  the  gaussian  aperture  moves  close  to  the  plane  mirror  (x  -*  0)  the  self- 
consistent  mode  approaches  a  spherical  wave  centered  on  the  plane  mirror, 
with  zero  spot  size  at  the  aperture  and  infinite  spot  size  at  the  PCM,  as 
in  Fig.  5(c) .  In  the  opposite  limit  as  the  gaussian  aperture  moves  close 
to  the  PCM  (X  -» l)  the  mode  becomes  essentially  a  half-confocal  resonator 
mode,  with  a  minimum  spot  size  w  =  (2Lx/tt)*^  at  the  PCM  end  as  in 
Fig.  5(d). 


B.  Plane-Mirror  Cavity  With  Gaussian  Aperture  and  Gaussian  Duct 

Consider  next  a  plane-mirror  cavity  of  length  L  having  a  gaussian  aperture 
located  immediately  in  front  of  (or  possibly  inside)  the  PCM,  but  with  the  re¬ 
mainder  of  the  cavity  also  filled  with  a  gaussian  duct,  or  a  transversely  vary¬ 
ing  gain  medium,  as  in  Fig.  6(a).  The  propagation  constant  k  in  this  medium 
is  assumed  to  have  a  complex  radial  variation  given  by 


k2(r)  =  k2  [1  -  jnr2] 


(2>) 
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The  case  >  0  corresponds  to  a  radially  Increasing  loss  In  the  medium,  or 
alternatively  to  a  gain  medium  with  a  maximum  gain  on  axis  and  radially  de¬ 
creasing  gain  off  axis.  Alternatively  ^  <  0  corresponds  to  an  "anti-duct" 
with  a  gain  minimum  (or  loss  maximum)  on  axis. 

If  we  use  the  same  definition  F  -  7Ta2/LX  as  above  for  the  gaussian  ap¬ 
erture,  then  the  matrix  elements  are 


A  «  1  -  2j  (l/F  +  tiL2) 
B  =  2L  (1  -  2jr]L2/3) 


(24) 


In  the  limit  of  a  weak  gaussian  aperture  (|f|  -*«>)  and  weak  gaussian  ducting 
(  |t}  |  -*  0)  the  self-consistent  solutions  are 


where  Y  is  the  parameter 

Y  —  Ft)L2  (26) 


As  J F I  -» «o  and  |  r|  |  -»  0  ,  the  parameter  Y  can  take  on  any  value  in  the 
range  -  w  £  Y  £  00  depending  on  the  relative  magnitudes  of  F  and  t| 
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That  Is  to  say*  Y  is  a  measure  of  the  relative  strengths  of  the  gaussian 
aperture  and  the  gaussian  duct  even  vhen  both  become  very  weak. 

A  confined  self-consistent  gaussian  mode  only  exists  for  values  of  Y 
within  the  ranges 


Y  =  FqL2  a  -  1  (27) 

or 

Y  =  Fiji2  £  -  3  (28) 

These  ranges  define  a  complete  semicircle  in  the  l/q  plane,  as  illustrated 
in  Fig.  6.  However  the  perturbation  stability  criterion  of  Eq.  (18)  is 
satisfied  only  for 


/  1  2nL2  \ 

(; +  — )  * 0  <s) 

For  F  >  0  ,  which  means  a  positive  or  converging  gaussian  aperture,  the 

overall  system  can  still  be  stable  even  if  the  gaussian  duct  becomes  an  "anti¬ 
duct"  with  radially  increasing  transmission,  but  this  holds  only  over  the 
limited  range  Y  =  FnL2  a  -  2/^  .  Alternatively  for  F  <  0  ,  which  would 

mean  a  negative  or  diverging  gaussian  aperture,  the  gaussian  duct  must  have 
a  transmission  maximum  on  axis  of  relative  strength  sufficient  to  make 
Y  *»  FqL  i  +  2/j5  .  Figure  6(b)  uses  dashed  lines  to  indicate  the  portions 
of  the  self-consistent  locus  that  represent  confined  but  not  perturbation- 
stable  solutions  for  the  case  F  >  0  .  This  example  rather  clearly 
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illustrates  how  a  mixture  of  positive  and  negative  gaussian  aperturing  can 
lead  to  such  confined  but  perturbation-unstable  modes. 


C.  General  Cavity  With  Gaussian  Aperture 


As  a  final  example  we  consider  a  PCM  resonator  with  a  gaussian  aperture 
immediately  in  front  of  the  PCM  and  an  arbitrary  real  ABCD  system  outside 
this.  The  complex  ABCD  values  then  become 


A 


B  <*  B 

r 


(30) 


where  Af  and  B^  are  the  purely  real  values.  The  self-consistent  solution 
then  becomes 


(31) 


For  an  arbitrarily  weak  gaussian  aperture  (X/rra^  -»  0)  the  self-consistent 
value  approaches  in  the  limit  the  point  at  the  top  of  the  semicircle  in  the 
B/q  plane.  No  matter  ho*?  weak  the  aperture,  the  mode  still  attempts  to 
minimize  the  spot  size  wgc  at  the  plane  of  the  aperture. 
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A  particular  case  of  interest  is  a  cavity  of  this  type  consisting  only  of 
a  mirror  of  radius  Rq  at  distance  L  from  the  PCM  and  the  gauss ian  aperture. 
The  real  matrix  elements  are  then 


Ar  «  1  -  2L/R0 

Br  *  2L( 1  -  l/rq) 


(32) 


Figure  7  illustrates  the  locus  of  self-consistent  solutions  for  different  values 
of  the  mirror  radius  Rq  in  the  limit  of  weak  aperturing,  plotted  here  in  the 
complex  L/q  plane,  not  the  complex  fi/q  plane.  All  of  these  solutions  are 
stable. 
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IV.  RESONANT  FREQUENCIES  OF  PCM  CAVITIES 


Our  final  objective  is  to  explore  the  resonant  frequency  behavior  and  the 
axial  modes  of  a  laser  cavity  having  a  PCM  reflector  at  one  end. 

A.  Frequency  Dependence  of  PCM  Reflectivity 

The  resonance  properties  of  an  optical  cavity  with  a  PCM  reflector  at  one 
end  will  depend  in  the  most  general  case  on  the  phase  versus  frequency  charac¬ 
teristics  of  the  PCM  reflector*  as  well  as  the  remainder  of  the  optical  cavity 
A  signal  wave  at  a  detuned  frequency  ci^  +  An  sent  into  a  PCM  consisting  of 
a  four-wave  mixing  cell  pumped  at  frequency  will  be  reflected  as  a  con¬ 

jugate  idler  wave  at  the  oppositely  detuned  frequency  u>q  -  An  .  In  addition* 
the  amplitude  of  the  reflection  coefficient  will  decrease  for  off-resonance 

signals  because  of  unavoidable  phase  mismatch  between  the  signal*  pump  and 

9 

idler  frequencies  over  the  finite  cell  length.  The  bandwidth  of  a  PCM 
reflector  by  itself*  in  the  weak-coupling  limit*  can  be  derived  as  follows. 

Consider  a  four-wave  mixing  cell  having  total  interaction  length  * 

with  the  reference  plane  taken  to  be  at  the  front  (i.e.  input)  plane  of  the 
cell.  Let  the  signal  input  be  an  off-resonance  unit  impulse*  that  is*  an 
optical  pulse  with  carrier  frequency  cn^  +  An  and  duration  much  shorter 
than  the  transit  time  l-c/c  through  the  cell.  As  this  pulse  travels 
through  the  four-wave  mixing  cell  in  the  +z  direction  and  interacts 
with  the  assumed  cw  pump  waves*  it  will  shed  a  continuous  "tail"  of 
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constant-amplitude  idler  waves  at  frequency  -  Au  traveling  back  in  the 
-a  direction  through  the  cell.  This  idler  radiation  will  in  fact  continue 
to  emerge  back  out  through  the  front  face  of  the  cell  for  one  full  transit 
time  after  the  signal  impulse  leaves  the  cell  at  the  opposite  end. 

The  impulse  response  of  the  four-wave  cell  is  thus  essentially  a  square 
pulse  of  constant  amplitude  and  time  duration  2L c/c  ,  extending  from  t  =  0 
to  t  «=  2L Jc  for  a  reference  plane  at  the  front  face  of  the  cell.  Shifting 
the  reference  plane  to  the  center  of  the  cell  makes  the  impulse  response 
symmetric  from  -  Lc/c  to  Lc/c  .  The  frequency  response  of  the  cell  is 
then  the  Fourier  transform  of  this  time  impulse  response;  and  the  power  re¬ 
flectivity  is  the  magnitude  of  this  frequency  response  squared,  or 


K(o>) 


sin  AuL  Jc 
_ c 

AoLc/c 


(33) 


This  agrees  with  the  results  derived  by  others'*  in  the  weak-coupling  limit. 
The  power  reflectivity  has  a  maximum  value  R(u^)  on  resonance  and  falls 
off  with  the  familiar  |(sin  x)/x|^  dependence  on  frequency  detuning  on 
either  side. 


B.  PCM  Cavity  Resonance  Condition 

The  PCM  reflector  thus  has  a  frequency-dependent  reflectivity  by  itself. 
However,  if  the  PCM  cell  length  is  small  compared  to  the  length  L  of 

the  remainder  of  the  optical  cavity,  then  the  central  lobe  of  the  PCM  fre¬ 
quency  response  will  be  broad  compared  to  the  optical  cavity’s  axial  mode 
spacing.  In  addition,  if  the  PCM  reference  plane  is  taken  at  the  center  of 
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the  four-wave  mixing  cell*  the  reflectivity  at  least  for  the  weak-coupling 
case  will  have  constant  absolute  phase  versus  frequency  across  the  main  fre¬ 
quency  lobe.  We  will  therefore  evaluate  the  resonance  properties  of  a  PCM 
cavity  assuming  a  flat  and  constant-phase  reflectivity  for  the  PCM  reflector. 
These  conditions  may  or  may  not  be  achieved  in  different  experimental  situations. 
However,  if  the  cell  length  is  not  short  compared  to  the  cavity  length  L  , 

then  the  off-resonance  axial  modes  will  be  so  far  down  on  the  reflectivity  curve 
of  the  PCM  that  they  may  have  little  chance  of  oscillating  in  competition  with 
the  central  mode  in  any  case.  This  latter  situation  has  been  analyzed  in  more 
detail  by  AuYeung  £t  al. 

The  round-trip  phase  shift  of  the  remainder  of  the  optical  cavity  will  also 
depend  slightly,  in  second  order,  on  the  transverse  mode  character  of  the  cavity 
because  of  the  well-known  Guoy  phase  shifts.  These  effects  cause  slight  trans¬ 
verse-mode  shifts  for  higher  order  transverse  modes,  which  we  are  not  con¬ 
sidering;  but  otherwise  produce  in  effect  only  a  minuscule  change  in  the 

V 

effect  length  of  the  cavity.  We  will  therefore  ignore  these  and  use  an 
essentially  transmission-line  or  plane-wave  approach  for  the  main  part  of 
the  optical  cavity. 

We  therefore  adopt  a  model  as  shown  in  Fig.  8,  consisting  of  a  four-wave 
mixing  cell  located  at  z  «  0  and  a  normal  mirror  located  at  z  «  -  L  along 
the  z  axis.  The  four-wave  mixing  cell  is  pumped  by  two  contradirected  waves 

IN/  /V 

of  complex  phasor  amplitude  E^  and  Eg  at  the  driving  frequency  u^ 

There  is  also  incident  on  the  cell  from  inside  the  cavity  a  signal  wave  of 
the  general  form 

^  Jov^t 

e(t)*=Rcn(t)c  ( _>!i ) 
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He  now  ignore  transverse  variations  (for  simplicity)  but  assume  that  both  the 
phase  and  amplitude  of  the  complex  phasor  E(t)  may  be  time  varying  in  general. 
We  also  assume  that  the  absolute  phases  of  all  waves  are  measured  at  a  ref¬ 
erence  plane  located  at  z  *»  0  ,  in  the  center  of  the  four-wave -mixing  cell. 

For  small  incident  fields,  this  will  produce  a  four-wave-ref lected  or  PCM  wave 
given  by 

/  v  .  .  H*  /  \ 

C'(t)  =  Re  KE  (t)  e  (35) 

The  four-wave-mixing  coefficient  K  (which  may  in  general  be  complex,  with  a 
fixed  phase  angle)  depends  on  the  cell  geometry  and  on  the  magnitude  and  phase 
of  the  nonlinear  susceptibility.  It  is  also  directly  proportional  to  the 

/v  j 

pumping  wave  phasor  amplitudes  E^  and  Eg  . 

If  the  round-trip  transit  time  around  the  cavity  from  the  PCM  to  the  normal 
mirror  and  back  is  T  =  2L/c  ,  the  waves  C(t)  and  G'(t)  will  also  be 
related  by 

C(t  +  T)  »  g1/2  e'(t)  (36) 

where  G  is  the  net  round-trip  power  gain  due  to  any  laser  material  (or 

losses)  inside  the  cavity.  For  simplicity  we  assume  G  has  no  dispersion. 

Since  we  are  interested  primarily  in  the  resonant  frequencies  of  the  cavity 

we  will  assume  that  this  round-trip  gain  and  the  PCM  reflectivity  combine  to 

i  l/s  — 1 

give  unity  net  gain,  i.e.  (G  '  K|  s  1  .  If  we  denote  the  net  phase  angle 


73  - 


of  the  quantity  k  by  ♦  ,  then  from  (35)  and  (36)  the  round-trip  self- 
consistency  condition  for  the  cavity  becomes 


joWT  Jifr  ^ 

E(t  +  T)  e  ^  «  e  E'(t) 


e  E  (t) 


Consider  first  the  case  when  E(t)  is  constant.  If  we  write  the  phasor  amp¬ 


litude  as 


E  =  Eq  e 


then  the  self-consistency  condition  becomes 


J(0O  +  V} 


■  Eo  e 


J(t  -  e0) 


The  cavity  has  a  resonance  at  cj^  independent  of  cavity  length  L  ,  but  the 
signal  at  this  frequency  must  have  a  definite  phase  0q  (measured  at  z  «  0) 


given  by 


3q  “  iff  -  o^T  (mod  2n) 


The  waves  are  in  essence  not  only  frequency-synchronized  but  also  phase- 

synchronized  to  the  frequency  of  the  pump  waves  as  modified  by  the 

«■>»* 

phase  angle  of  K 

More  generally,  suppose  that  c(t)  contains  an  upshifted  signal  component 

at  some  frequency  oc^  +  uv  where  <nm  is  some  as  yet  unspecified  modulation 

frequency.  The  phase  of  E(t)  then  varies  as 

Ju>  t 

E(t)«ce  m  (111) 
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This  will  produce  a  reflected  wave 


E'(t)  ec  E*(t) 


”Jo>  t 

in 

«  e 


(42) 


An  upshifted  component  of  £(t)  at  thus  produces  a  downshifted 

component  of  C'(t)  at  j  which  will  eventually  travel  around  to 

become  £(t)  .  We  must  therefore  write  E(t)  as  the  combination 


E(t) 


jcD  t 
~  m 

c+  e  + 


-ju)  t 
m 


(43) 


where  the  c's  are  complex  phasor  amplitudes.  Putting  this  into  the  fre¬ 
quency  self-consistency  relation  (3  )  and  separating  like- frequency  terms 
gives 


c+  =  e 


-ju>  T  „ 
m 


c+ 


(44) 


Combining  the  first  of  these  relations  with  the  complex  conjugate  of  the  second 
then  produces  the  condition 

j&o  T  jn2rr 

01 


where  n  is  any  integer.  The  modulation  frequency  or  sideband  frequency  cnm 
can  have  any  of  the  values  given  by 


Those  are  spaced  by  half  of  the  usual  axial  mode  spacing  for  a  cavity  of 
length  L  • 

The  field  C(t)  In  the  cavity  may  thus  be  expanded  in  resonant  modes 
In  the  form 


e(t) 


E(t) 


j(o)gt  +0q) 


ZS 


jn(7t/T)t 
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j^t  +0q) 
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with  the  condition  that 


c  =  (-1)  c 

-n  '  n 


(48) 


This  is  equivalent  to  writing 


E(t)  *=  Eg  +  j  En  sin  (n7Tt/x  +  0n) 


n=l,2,  . 


^  En  cos  (irnt/r  +  0n) 

n-2, 4 . . . 


(49) 


The  PCM  cavity  has  a  central  resonance  frequency  determined  by  the  pumping 
frequency  (not  the  cavity  length  h) .  This  central  frequency  is  sur¬ 
rounded  by  a  set  of  paired  "half-axial-modc "  resonances  spaced  by  c/4l  llz  , 
or  half  the  usual  spacing  for  a  cavity  of  length  L 
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This  behavior  may  be  understood  physically  by  the  following  example. 
Assume  that  at  an  initial  instant  the  cavity  is  entirely  filled  with 
radiation  at  an  unshifted  frequency  +  M  ,  as  in  Fig.  9(a)  which  shows 
the  waves  traveling  in  both  directions  inside  the  cavity.  The  radiation 
moving  toward  the  PCM  will  be  reflected  at  the  downshifted  frequency 
uiq  -  An  ,  and  therefore  after  a  short  time  the  signal  inside  the  cavity 
will  appear  as  in  Fig.  9(b),  with  lower-frequency  radiation  now  moving 
away  from  the  PCM.  The  PCM  will  continue  to  take  in  the  higher  frequency 
o)q  +  Ad  and  emit  the  lower  frequency  o>q  -  An  until  after  one  complete 
round-trip  time  T  the  cavity  will  be  entirely  filled  with  lower  frequency 
radiation.  But,  the  process  will  then  reverse  and  the  incident  signal  at 
u)q  -  An  will  be  converted  back  to  the  original  value  +  An  •  This 
cycle  will  be  complete  after  a  time  2T  ,  or  two  round  trips.  It  is  the 
doubling  of  the  time  interval  before  the  signal  returns  to  its  original 
state  that  leads  to  the  halving  of  the  effective  axial  mode  spacing. 

Although  the  apparent  or  instantaneous  frequency  of  the  cavity  signal 
can  oscillate  back  and  forth  about  the  central  value  mg  with  period  2T  , 
amplitude  variations  will  still  self-reproduce  after  period  T  .  Mode 
locking  of  this  cavity,  for  example,  will  produce  pulses  at  intervals  of 
T  ,  not  2T  .  However,  the  apparent  frequency  of  every  alternate  pulse 
could  oscillate  back  and  forth  about  cdq  with  period  2T 


V.  CONCLUSIONS 


We  have  carried  out  the  gaussian  transverse  mode  analysis  of  a  resonator 
formed  by  conventional  paraxial  optical  components  plus  a  phase  conjugate 
mirror  (PCM).  For  a  resonator  with  only  real  paraxial  elements  (no  transverse 
amplitude  filtering)  any  gaussian  beam  reproduces  itself  after  two  round  trips 
and  there  is  a  continuous  set  of  self-consistent  q  values  that  reproduce 
themselves  after  one  round  trip.  These  q  values  correspond  to  all  possible 
gaussian  beams  which  match  the  mirror  curvature  but  have  arbitrary  spot  size> 
at  the  opposite  end  of  the  cavity  from  the  PCM. 

If  complex  paraxial  elements  such  as  gaussian  apertures  are  added  to  the 
resonator*  there  is  only  one  stable  self-reproducing  or  self-consistent 
gaussian  mode  of  the  resonator.  Gaussian  beams  starting  from  any  other  point 
in  the  l/q  plane  converge  to  this  self-consistent  value  upon  repeated  round 
trips . 

The  resonant  frequency  structure  of  the  PCM  cavity  contains  a  central 
resonance  exactly  at  the  PCM  pumping  frequency  *  independent  of  the 

cavity  length  L  .  This  frequency  is  surrounded  by  a  set  of  "half-axial" 
resonant  modes  spaced  by  c/^L  ,  or  half  the  usual  c/2L  axial  mode  spacing 
for  a  cavity  of  length  L  .  However*  there  are  phase  and  amplitude  con¬ 
straints  on  the  phasor  amplitudes  of  each  pair  of  sidebands  above  and  below 
the  central  frequency. 

This  work  was  supported  by  the  U.S.  Air  Force  Office  of  Scientific 
Ho search . 
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CAPTIONS  FOR  FIGURES 


1.  (a)  A  typical  PCM  resonator,  and  (b)  the  basic  analytical  model  for 
discussion  of  the  lowest-order  gaussian  transverse  mode  of  an  optical 
cavity  with  a  phase-conjugate  mirror.  The  gaussian  q  parameters  are 
measured  at  a  reference  plane  immediately  in  front  of  the  PCM. 

2.  (a)  Locus  of  self-consistent  gaussian  beam  modes  in  the  |B | /q 
plane  for  real  ABCD  elements,  (b)  A  single  round-trip  inside  the 
PCM  cavity  transforms  an  arbitrary  inital  point  B/q  to  the  geometri¬ 
cally  inverted  point  B/q'  with  respect  to  the  displaced  unit  circle. 
A  second  round  trip  returns  B/q'  to  B/q 

Unfolded  model  for  two  successive  round  trips  in  a  simple  PCM  cavity. 
The  PCM  acts  like  a  converging  lens  for  a  diverging  wave  (solid  lines) 
but  like  a  diverging  lens  for  a  converging  wave  (dashed  lines). 

4.  A  PCM  cavity  with  weak  gaussian  aperturing  has  a  unique  lowest-order 
self-consistent  gaussian  mode  with  a  gaussian  parameter  qgc  ,  cor¬ 
responding  to  a  single  point  located  on  or  close  to  the  displaced 
unit  circle.  A  beam  starting  from  any  other  initial  value  of  q 

converges  toward  the  value  q  on  repeated  bounces . 

s  c 

5.  Example  of  a  complex  PCM  optical  cavity  containing  a  plane  mirror 
(PM),  gaussian  aperture  (GA),  and  phase-conjugate  mirror  (PCM) . 
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(a)  Cavity  model,  (b)  Self-consistent  points  in  the  l/q  plane 

for  X  «  L  /L  ranging  from  0  to  1  .  (c)  Self-consistent  mode  for 

X  -> 0  ;  the  spot  size  approaches  zero  at  one  end  and  infinity  at  the 
other,  (d)  Self-consistent  mode  for  X  -»  1  ;  the  mode  is  essentially 
a  half-confocal  resonator  mode. 

6.  (a)  Optical  cavity  containing  a  plane  mirror  (PM),  gaussian  duct  (GD), 

gaussian  aperture  (GA),  and  phase-conjugate  mirror  (PCM) .  (b)  Locus 

of  self-consistent  solutions  in  the  l/q  plane  for  various  values 

2 

of  the  relative-strength  parameter  Y  =  Fi)L  .  Solutions  on  the  solid 
line  are  perturbation-stable  while  those  along  the  dashed  regions  are 
confined  but  perturbation-unstable  (for  the  case  F  >  0). 

7*  Locus  of  stable  self-consistent  solutions  for  a  cavity  length  L 

containing  a  curved  mirror  of  radius  Rq  ,  a  weak  gaussian  aperture, 
and  a  phase  conjugate  mirror,  with  different  values  of  the  parameter 

r0/l  . 

8.  Model  for  analysing  the  resonant  frequencies  of  a  PCM  cavity  of  length 
L  .  The  PCM  cell  is  pumped  at  frequency  cdq  and  has  a  reference 
plane  at  z  =  0  . 

9*  (a)  A  PCM  cavity  is  initially  filled  with  radiation  at  u>q  +  An 

(b)  Reflection  from  the  PCM  converts  the  incident  radiation  at 
cuq  +  An  to  coq  -  An  ,  as  shown  here  after  ~  l/4  of  a  cavity 
round  trip,  (c)  After  ~  1  complete  round  trip,  the  cavity  is 
now  nearly  filled  with  radiation  at  tug  -  An  , and  the  reverse 
process  begins . 
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FIGURE  2(b) 
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